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Abstract 

Punctured polygons are polygons with internal holes which are also polygons. 
The external and internal polygons are of the same type, and they are mutually as 
well as self-avoiding. Based on an assumption about the limiting area distribution 
for unpunctured polygons, we rigorously analyse the effect of a finite number of 
punctures on the limiting area distribution in a uniform ensemble, where punctured 
polygons with equal perimeter have the same probability of occurrence. Our analysis 
leads to conjectures about the scaling behaviour of the models. 

We also analyse exact enumeration data. For staircase polygons with punctures 
of fixed size, this yields explicit expressions for the generating functions of the first 
few area moments. For staircase polygons with punctures of arbitrary size, a careful 
numerical analysis yields very accurate estimates for the area moments. Interest- 
ingly, we find that the leading correction term for each area moment is proportional 
to the corresponding area moment with one less puncture. We finally analyse cor- 
responding quantities for punctured self-avoiding polygons and find agreement with 
the conjectured formulas to at least 3-4 significant digits. 

1 Introduction 

The behaviour of planar self-avoiding walks (SAW) and polygons (SAP) is one of the clas- 
sical unsolved problems, not only of algebraic combinatorics, but also of chemistry and of 
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physics [H [21 [3] . In the field of algebraic combinatorics, it is a classical enumeration prob- 
lem. In chemistry and physics, SAWs and SAPs are used to model a variety of phenomena, 
including the properties of long-chain polymers in dilute solution [4], the behaviour of ring 
polymers and vesicles in general [5] and benzenoid systems P, [7] in particular. Though the 
qualitative form of the phase diagram [8] is known rigorously, there is otherwise a paucity 
of rigorous results. However, there are a few conjectures, including the exact values of 
the critical exponents [9l[T0], and more recently the limit distribution of area and scaling 
function for SAPs, when enumerated by both area and perimeter [Til 121 EES EEH [15] - 

Models of planar polygons with punctures arise naturally as cross-sections of three- 
dimensional vesicle models. In such cross-sections, there may be holes within holes, and 
the number of punctures may be infinite. In this work, we exclude these possibilities. 
Whereas our methods can be used to study the former case, the second situation presents 
new difficulties, which we have not yet overcom^j]. 

In this work we consider the effect of a finite number of punctures in polygon models, 
in particular we study staircase polygons and self-avoiding polygons on the square lattice. 
The perimeter of a punctured polygon [161 E] is the perimeter of its boundary (both 
internal and external) while the area of a puncturedpolygon is the area of the enclosed by 
the external perimeter minus the area(s) of any holesj. As discussed in section [2] below, the 
effect of punctures on the critical point and critical exponents of the area and perimeter 
generating function has been the subject of previous studies, but the effect of punctures 
on the critical amplitudes and detailed asymptotics have not, to our knowledge, been 
previously considered. 

Apart from the intrinsic interest of the problem, we also believe it to be the appropriate 
route to study the detailed asymptotics of polyominoes, since punctured polygons are a 
subclass of polyominoes. While we still have some way to go to understand the polyomino 
phase diagram, we feel that restricting the problem to this important subclass is the correct 
route. 

The make-up of the paper is as follows: In the next section we review the known 
situation for the perimeter and area generating functions of punctured polygons and poly- 
ominoes. In section 3 we review the phase diagram and scaling behaviour of staircase 
polygons and self-avoiding polygons. In section 4 we rigorously express the asymptotic 
behaviour of models of punctured polygons in the limit of large perimeter in terms of 
the asymptotic behaviour of the model without punctures, by refining arguments used in 
|16j . This leads, in particular, to a characterisation of the limit distribution of the area of 
punctured polygons. This result is then used to conjecture scaling functions of punctured 

1 Since punctured polygons with an unlimited number of punctures have, in contrast to polygons without 
punctures, an (ordinary) perimeter generating function with zero radius of convergence 18J, both the phase 
diagram and the detailed asymptotics are clearly going to be very different from those of polygons without 
punctures. This is discussed further in the conclusion. 

2 This has to be distinguished from so-called composite polygons [IS]. The perimeter of a composite 
polygon is defined as the perimeter of the external polygon only, resulting in asymptotic behaviour different 
from punctured polygons. Moreover, composite polygons can have more complex internal structure than 
just other polygons. 
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polygons. We consider three cases of increasing generality. First, we consider the case of 
minimal punctures. It is shown that effects of self- avoidance are asymptotically irrelevant, 
and that elementary area counting arguments yield the leading asymptotic behaviour. We 
then discuss the case of a finite number of punctures of bounded size, and finally the case 
of a finite number of punctures of unbounded size. Results for the latter case are given 
for models with a finite critical perimeter generating function such as staircase polygons 
and self-avoiding polygons. Whereas the latter two cases are technically more involved, the 
underlying arguments are similar to the case of minimal punctures. If the critical perimeter 
generating function of the polygon model without punctures is finite, then all three cases 
lead, up to normalisation, to the same limit distributions and scaling function conjectures. 

The next two sections discuss the development and application of extensive numerical 
data to test the results of the previous section. Moreover, the numerical analysis yields 
predictions, conjectured to be exact, for the corrections to the asymptotic behaviour. In 
particular, section 5 describes the very efficient algorithms used to generate the data, while 
section 6 applies a range of numerical tools to the analysis of the generating functions for 
punctured staircase polygons and then punctured self-avoiding polygons. Here we wish to 
emphasise that our work on this problem involved a close interplay between analytical and 
numerical work. Initially, our intention was to check our predictions for scaling functions 
by studying amplitude ratios for area moments (given in Tabled]). We subsequently dis- 
covered numerically the exact solutions for minimally punctured staircase polygons. We 
also obtained very accurate estimates for the amplitudes of staircase polygons with one 
or two punctures of arbitrary size. From these results we were able to conjecture exact 
expressions for the amplitudes, which in turn spurred us on to further analytical work in 
order to prove these results. The final section summarises and discusses our results. 

2 Punctured polygons 

We consider polygons on the square lattice in this article. In particular, we study self- 
avoiding polygons and staircase polygons. A self-avoiding polygon on a lattice can be 
defined as a walk along the edges of the lattice, which starts and ends at the same lattice 
point, but has no other self- intersections. When counting SAPs, they are generally con- 
sidered distinct up to translations, change of starting point, and orientation of the walk, 
so if there are p m SAPs of length or perimeter m there are 2mp m walks (the factor of two 
arising since the walk can go in two directions). On the square lattice the perimeter of 
any polygon is always even so it is natural to count polygons by half-perimeter instead 
of perimeter. The area of a polygon is the number of lattice cells (times the area of the 
unit cell) enclosed by the perimeter of the polygon. A (square lattice) staircase polygon 
can be defined as the intersection of two mutually avoiding directed walks starting at the 
same lattice point, moving only to the right or up and terminating once the walks join 
at a vertex. Every staircase polygon is a self-avoiding polygon. It is well known that the 
number p m of staircase polygons of half-perimeter m is given by the (m — l) th Catalan 
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Figure 1: Examples of the types of staircase polygons we consider in this paper, 
number, p m = P"^ 2 ) /m, with half-perimeter generating function 



V(x) = 



PrnX 
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where the connective constant \i — 4 and the critical exponent a = 3/2. Recall that 
f(x) ~ g(x) as x S x c means that \im f(x)/g(x) = 1 as x — > x c from below. In addition, 
as usual, the rhs is understood as the first two leading terms in an asymptotic expansion 
of the lhs about x = l//i, see e.g. [201 Sec 1]. 

Punctured polygons [IB] are polygons with internal holes which are also polygons (the 
polygons are mutually- as well as self-avoiding). The perimeter of a punctured polygon is 
the sum of the external and internal perimeters while the area is the area of the external 
polygon minus the areas of the internal polygons. We also consider polygons with minimal 
punctures, that is, polygons where the punctures are unit cells (or polygons with perimeter 
4 and area 1). Punctured staircase polygons are illustrated in figure [U 

We briefly review the situation for SAPs with punctures. Analogous results can be 
shown to hold for staircase polygons with punctures. Square lattice SAPs with r punctures, 
counted by area n, were first studied by Janse van Rensburg and Whittington [17]. They 
proved the existence of an exponential growth constant satisfying 
Denoting the corresponding number of SAPs by On and assuming asymptotic behaviour 
of the form 



^(r)( K M)n n/ 3 r -l ( 



n 



oo 



Janse van Rensburg proved [21] that (5 r = j3o + r. These results of course translate to 
the singular behaviour of the corresponding generating functions, defined by A^ r \q) = 

In [16J Guttmann, Jensen, Wong and Enting studied square lattice SAPs with r punc- 
tures counted by half-perimeter m. They proved the existence of an exponential growth 
constant /i^ satisfying jj,^ = /i^ -* = /x. If the corresponding number of SAPs is 



assumed to behave asymptotically as 
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they argued, on the basis of a non-rigorous argument, that a r = «o + |r. Their results 

also translate to the associated half-perimeter generating function V^(x) = ^2 m>0 PmX m 
corresp ondingly. 

Similar results were obtained for polyominoes enumerated by number of cells (i.e. area) 
with a finite number r of punctures [16]. It has been proved that an exponential growth 
constant r exists independently of r, which satisfies 4.06258 ~ r > k « 3.97087, where k 
is the growth constant for SAPs enumerated by area. If the number a„ of polynominoes 
of area n with r punctures is assumed to satisfy asymptotically 

a (r) _ ^{t^V 1 (n -> oo), 

it has been shown that 7 r = 7o + r and hence that, if the exponents j r exist, they increase 
by 1 per puncture. It was further conjectured on the basis of extensive numerical studies 
[15] . that the number an satisfies asymptotically 

<4 r) ~ r n n r ~ l C^/n* (n -> oo). 

i>0 

Notice the conjecture 70 = and that the correction terms go down by a whole power. 

For unrestricted polyominoes, that is to say, with no restriction on the number of 
punctures, it was proved by Guttmann, Jensen and Owczarek [18J that the perimeter 
generating function has zero radius of convergence. The perimeter is defined to be the 
perimeter of the boundary plus the total perimeter of any holes. If p m denotes the number 
of polyominoes, distinct up to a translation, with half-perimeter m, they proved that 
p m = m m / 4 +°( m ) ; meaning that 

\ogp m 1 
hm — = -. 

nwoo m log m 4 

An attempt to study the quasi-exponential generating function with coefficients r m = 
p m /T(m/A + 1) was equivocal. For that reason, studying punctured self-avoiding polygons 
was considered a controlled route to attempt to determine the two-variable area-perimeter 
generating function of polyominoes. 

In passing, we note that in [22] the exact solution of the perimeter generating function 
for staircase polygons with a staircase hole is conjectured, in the form of an 8 th order 
ODE. It is not obvious how to extract particular asymptotic information, notably critical 
amplitudes from the solution without numerically integrating the ODE. In the following, 
we will obtain such information by combinatorial arguments, which refine those of [16J. 

3 Polygon models and their scaling behaviour 

We review the asymptotic behaviour of self-avoiding polygons and staircase polygons fol- 
lowing mainly [8 J . For concreteness, consider the fixed perimeter ensemble where, for fixed 
half-perimeter m, each polygon of area n has a weight proportional to q n , for some positive 
real number q. If < q < 1, polygons of large area are exponentially suppressed, so that 
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typical polygons should be ramified objects. Since such polygons would closely resemble 
branched polymers, the phase < q < 1 is also referred to as the branched polymer phase. 
As q approaches unity, typical polygons should fill out more, and become less string-like. 
For q > 1, polygons of small area are exponentially suppressed, so that typical polygons 
should become "fat" . Indeed, they resemble convex polygons [23] and it has been proved 
[8] that the mean area of polygons of half-perimeter m grows asymptotically proportional 
to m 2 . In the extended phase q — 1, it is numerically very well established that the mean 
area of polygons of half-perimeter m grows asymptotically proportionally to m 3 / 2 . In the 
branched polymer phase < q < 1, the mean area of polygons of half-perimeter m is ex- 
pected to grow asymptotically linearly in m, compare also [2H Thm 7.6] and [251 Ch IX. 6, 
Ex. 12]. 

This change of asymptotic behaviour of typical polygons w.r.t. q is reflected in the 
singular behaviour of the half-perimeter and area generating function 

V(x,q)=J2Pm,nX m q n , 

where p m ,n denotes the number of (self-avoiding) polygons of half-perimeter m and area n. 
It has been proved [8] that the free energy 



n(q) :-- 



hm — log [ 

\ n 



Pm,nq 



exists and is finite if < q < 1. Further, n(q) is log-convex and continuous for these 
values of q. It is infinite for q > 1. It was proved that for fixed < q < 1, the radius of 
convergence x c (q) of V(x,q) is given by x c (q) = e~ K<yq \ For fixed q > 1, V(x,q) has zero 
radius of convergence. Fisher et al. [8] obtained rigorous upper and lower bounds on x c (q). 
The expected phase diagram, i.e., the radius of convergence of V(x, q) in the x — q plane, 
as estimated numerically from extrapolation of SAP enumeration data by perimeter and 
area, is sketched qualitatively in figure [2] 




1 q 



Figure 2: A sketch of the phase diagram of self- avoiding polygons. 
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For < q < 1, the line x c (q) is, for self-avoiding polygons, expected to be a line of 
logarithmic singularities of the generating function V(x, q). For branched polymers in the 
continuum limit, the existence of the logarithmic singularity has recently been proved [26] . 
The line q — 1 is, for < x < x c := x c (l), a line of finite essential singularities [8J. For 
staircase polygons, counted by half-perimeter and area, the corresponding phase diagram 
can be determined exactly, and is qualitatively similar to that of self-avoiding polygons. 
Along the line x c (q) the half-perimeter and area generating function diverges with a simple 
pole, and the line q — 1 is, for < x < x c , a line of finite essential singularities 



We will focus on the uniform fixed perimeter ensemble q — 1 in this article. Whereas 
asymptotic area laws in the fixed perimeter ensemble are expected to be Gaussian for posi- 
tive q 1, the behaviour in the uniform fixed perimeter ensemble q — 1 is more interesting. 
For staircase polygons, it can be shown that a limit distribution of area exists and is given 
by the Airy distribution [28], [29], [30]. For self-avoiding polygons, it is conjectured that 
an area limit law exists and is given by the Airy distribution, on the basis of a detailed 
numerical analysis [TTJ [HI [T5] . See subsections 14.11 and 14.41 

If Pm,n denotes the number of polygons of half-perimeter m and area n, the existence 
and the form of a limit distribution can be inferred from the asymptotic behaviour of the 
factorial moment coefficients J2 n {n)kPm,n, where (a)k = a ■ (a — 1) • . . . • (a — k + 1). The 
following result is obtained by standard reasoning |31j . 

Proposition 1. Let for m,n e No real numbers p m ,n be given. Assume that the numbers 
p m ,n have the asymptotic form, for k 6 Kf , 



n)kPm,n ~ A k x c m m rk 1 (m — ► oo) (2) 



for positive real numbers andx c , where 7& = (k — 6)/<f), with real constants 9 and (p > 0. 
Assume that the numbers := A^/Aq satisfy the Carleman condition 

oo 

Y^{M 2k )- l ' 2k = +oo. (3) 

k=0 



Then, for almost all m, the random variables X m of area in the uniform fixed perimeter 
ensemble 

= n) 



Pm,n 



Pm,r 



are well defined. We have 



X m := — — —> X (m — > oo) 



m 



1/4 



for a uniquely defined random variable X with moments M k , where the superscript d denotes 
convergence in distribution. We also have moment convergence. 
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Sketch of proof. A straightforward calculation using Eq. (|2J) leads to 

E[(X m ) k ] ~ ^™ fc/ * (™ - °°)- 

It follows that asymptotically the factorial moments are equal to the (ordinary) moments. 
Thus, the moments of X m have the same asymptotic form 

E[(X m ) k ] ~ ^ = M k (m-oo). 

Due to the growth condition Eq. Q, the sequence (M fc ) fcgN() defines a unique random 
variable X with moments M k . Moment convergence of (X m ) implies convergence in dis- 
tribution, see [311 Thm 4.5.5] for the line of arguments. □ 

The assumption Eq. (j2J) translates, on the level of the half-perimeter and area generat- 
ing function V(x,q), to a certain asymptotic behaviour of the so-called factorial moment 
generating functions 

( X \k Qk 

It can be shown (compare [32]) that the asymptotic behaviour Eq. (T5]) implies for 7*. > 
the asymptotic equivalence 

~ (Xc [ k x ^ k ( x y o> ( 4 ) 

where the amplitudes fk are related to the amplitudes in Proposition [1] via 

fk = { -^A k x^T( lk ). (5) 

If — 1 < 7fc < 0, the series g k (x) is convergent as 1 / x c , and the same estimate Eq. (J4]) 
holds, with ^(x) replaced by g k (x) — gk( x c), where g(x c ) := \im x y Xc g(x). In order to deal 
with these two different cases, we define for a power series g(x) with radius of convergence 
x c , the number 

( C ) = f g(x c ) if I lini^^x)! < 00 
1 otherwise. 

Adopting the generating function point of view, the amplitudes f k determine the numbers 
A k and hence the moments M k = A k /Ao of the limit distribution. The formal series 
F( s ) = J2k>o fk s ~ lk will appear frequently in the sequel. 



3 Note that our definition of the amplitudes A k differs from that in (T3] by a factor of (— l) k k\ and from 
that in [33l [34] by a factor of k\. 
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Definition 1. For the generating function V(x,q) of a class of self- avoiding polygons, 
denote its factorial moment generating functions by 



(-l) k d k 

Assume that the factorial moment generating functions satisfy 

9k(x) ~ gf ~ 7 ~ \Z7 (x / x c ), (6) 

with real exponents 7^. Then, the formal series 

F(s) = J2 hs~ lk 

fc>0 

is called the area amplitude series. 

The area amplitude series is expected to approximate the half-perimeter and area gen- 
erating function V(x, q) about (x, q) = (x c , 1). This is motivated by the following heuristic 
argument. Assume that 7^ = (k — 6) /(ft with (ft > and argue 

H*,*) « E fc > (si c) + (i - i) k 

- (E fc > ^ c) (i - a) k ) + (i - a) e {T. k >_oh (tw)^) ■ 

In the above calculation, we formally expanded V(x, q) about q — 1 and then replaced 
the Taylor coefficients by their leading singular behaviour about x = x c . In the rhs of the 
above expression, the first sum is by assumption finite, and the second term contains the 
area amplitude series F(s) of combined argument s = (x c — x)/(l — qy. This motivates the 
following definition. A class of self-avoiding polygons is a subset of self-avoiding polygons. 
Prominent examples are, among others [35J, self-avoiding polygons and staircase polygons. 

Definition 2. Let a class of square lattice self-avoiding polygons be given, with half- 
perimeter and area generating function V(x,q). Let < x c < 00 be the radius of conver- 
gence of the half-perimeter generating function V(x, 1). Assume that there exist a constant 
sq G [—00, 0), a function T : (s , 00) — > R, a real constant A and real numbers 9 and (ft > 0, 
such that the generating function V(x, q) satisfies, for real x and q, where < q < 1 and 
(x c — x)/{l — q)^ G (sqj 00), the asymptotic equivalence 

V(x, q)~A+{l- q) e T (tjt^) (*, ?) — (*« I)- (7) 

Then, the function J 7 (s) is called a scaling function of combined argument s = (x c — x)/(l — 
qy, and 9 and (ft are called critical exponents. 
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Remarks, i) Due to the restriction on the argument of the scaling function, the limit 
(x, q) — > (x c , 1) is approached for values (x, q) satisfying x < x (q) and q < 1, where 
Xo{q) = x c - s (l - q)*. 

ii) The above scaling form is also suggested by the theory of tri critical scaling, adapted 
to polygon models [36]. The scaling function describes the leading singular behaviour of 
V(x, q) about the point (x c , 1) where the two lines of qualitatively different singularities 
meet. 

Hi) The additional condition (f> > 9 and 6 ^ N ensure that 7^ G (— 1, 00) \ {0}. Then, 
by the above argument, it is plausible that there exists an asymptotic expansion of the 
scaling function J-'(s) about infinity coinciding with the area amplitude series F(s), i.e., 
JF(s) ~ F(s) as s — > 00. Recall that s is considered to be a real parameter. 

For staircase polygons the existence of a scaling form Eq. ([7]) has been proved [271 
Thm 5.3], with scaling function J-'(s) : (sq, 00) — > R explicitly given by 

^) = ^logAi(2 8 / 3 S ), (8) 

with exponents 9 = 1/3 and (p = 2/3 and x c = 1/4, where Ai(x) = - J °° cos(t 3 /3 + tx)dt is 
the Airy function. The constant so is such that 2 8 ^ 3 so is the location of the Airy function 
zero of smallest modulus. For rooted SAPs with half-perimeter and area generating function 
V r (x, q) = x-^T'(x, q), the conjectured form of the scaling function J>(s) : (so, 00) — > R is 

m 

^w = ^ 1oeA 'U( 4 ^ ss J- 

with the same exponents as for staircase polygons, 6 = 1/3 and <fi = 2/3. Here, x c = 
0.14368062927(2) is the radius of convergence of the half-perimeter generating function 
V r (x, 1) of (rooted) SAPs, and A = 0.09940174(4) is the critical amplitude ^2 n mPm,n ~ 
A x~ m m~ 3 / 2 of rooted SAPs, which coincides with the critical amplitude A of (unrooted) 
SAPs. Again, the constant Sq is such that the corresponding Airy function argument is 
the location of the Airy function zero of smallest modulus. This conjecture was based 
on the conjecture that both models have, up to normalisation constants, the same area 
amplitude series. The latter conjecture is supported numerically to very high accuracy 
by an extrapolation of the moment series using exact enumeration data [HI Q3]. The 
conjectured form of the scaling function J-"(s) : (so, 00) — > R for SAPs is obtained by 
integration, 

f(s) = - 7 ^\ogAi(-(4A )h)+C(q), (9) 
2tt \x c J 

with exponents 6 = 1 and = 2/3. In the above formula, C(q) is a q dependent constant of 
integration, C(q) = Tg-(1 — q) log(l — q), see [15]. Corresponding results for the triangular 
and hexagonal lattices can be found in [llj . 

For models of punctured polygons with a finite number of punctures, we have quali- 
tatively the same phase diagram as for polygon models without punctures, however with 
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different critical exponents 9 depending on the number of punctures [211 [16] , and hence we 
expect different scaling functions. We will focus on critical exponents and area limit laws 
in the uniform ensemble q — 1 in the following section. This will lead to conjectures for 
the corresponding scaling functions. 



4 Scaling behaviour of punctured polygons 

We briefly preview the main results of this section. In subsection 14.11 we study polygons 
with a finite number of minimal punctures. Our result assumes a certain asymptotic form 
for the area moment coefficients for unpunctured polygons. This 'assumed' form is known 
to be true for staircase polygons and many other models and universally accepted as true 
for self-avoiding polygons. Given this assumption, we prove that the asymptotic behaviour 
of the area moment coefficients for minimally punctured polygons can be expressed in terms 
of the asymptotic behaviour of unpunctured polygons. In particular we derive expressions 
for the leading amplitude of the area moments for punctured polygons in terms of the 
amplitudes for unpunctured polygons. For staircase polygons this leads to exact formulas 
for the amplitudes. For self-avoiding polygons the formulas contain certain constants which 
aren't known exactly but can be estimated numerically to a very high degree of accuracy. 
In subsection 14.21 we extend the study and proofs to polygons with a finite number of 
punctures of bounded size and then in subsection 14.31 to models with punctures of arbitrary 
or unbounded size. Finally in subsection 14.41 we consider the consequences of our results 
for the area limit laws of punctured polygons and we present conjectures for the scaling 
functions. 



4.1 Polygons with r minimal punctures 

For polygon models with rational perimeter generating functions, corresponding models 
with minimal punctures have been studied in [37] . In particular, a method to derive explicit 
expressions for generating functions of exactly solvable models with a minimal puncture 
was given [37J, Appendix]. It has been applied to Ferrers diagrams, whose perimeter and 
area generating function satisfies a linear g-difference equation, see J37, Eq. (54)]. The 
method can also be applied to the model of staircase polygons, whose half-perimeter and 
area generating function V(x, q) satisfies the quadratic g-difference equation 

^ ' ^ 1 — 2qx — V(qx, q) ^ 

Let V D ^ (x, q) denote the half-perimeter and area generating function of staircase polygons 
with r minimal punctures. We have the following result for the case r = 1. 

Fact 1. The half-perimeter and area generating function of staircase polygons with a single 
minimal puncture V D ^(x , q) is given by 

x A ( dV dV \ 

PD(1)(X ' g) = (l-2qx-V(qx,qW ^ q) - qx — (qx, q) + q — (qx, q) j , (11) 
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where V(x, q) satisfies Eq. /[TU\) . 



□ 



Remarks, i) For a proof of FactdJ proceed along the lines of [37J Appendix]. We do not 
give the details, since we are mainly interested in asymptotic results, for which we will give 
an elementary combinatorial derivation, valid for arbitrary r. See Proposition [2] and its 
subsequent extensions. 

ii) For polygons with r punctures, their k th area moment generating functions are defined 

by ^ (x) = (y~§^j *P ° ^ ( x 5 <?) I q -y The above equations can be used to obtain explicit 

expressions for the area moment generating functions V° ^ (x) by implicit differentiation. 

The functions \x) also appear in section ICTTl 

Hi) Assuming that V D ^'(x,q) has scaling behaviour of the form 



about (x, q) = (x c , 1), and the necessary analyticity conditions for the validity of the 
following calculation, we can express the scaling function J rD( - 1 \s) of staircase polygons 
with a single minimal puncture in terms of the known scaling function T{s) of staircase 
polygons Eq. ([8]). From Eq. (fTTj) we infer that 61 = —2/3, 0i = 2/3 and 



In principle, the method of [37J Appendix] can be used to analyse the case of several 
minimal punctures. However, the analysis becomes quite cumbersome. On the other 
hand, the previous result suggests simple expressions for the scaling functions of models 
with several punctures in terms of that without a puncture. Moreover, we expect such 
a phenomenon also to occur for models where an exact solution does not exist or is not 
known. This is discussed next. We will asymptotically analyse the area moments of a 
polygon model with punctures and draw conclusions about their possible scaling behaviour. 

For a class of punctured self-avoiding polygons, consider their area moment coefficients 



where p m ,n denotes the number of polygons in the class with r minimal punctures, r G No, 
of half-perimeter m and area n. For simplicity of notation, we write p mjT1 := pm,« and 
Pm '■= Pm°' k ^ ■ The area moments in the uniform fixed perimeter ensemble are expressed 
in terms of the area moment coefficients via 



V 



q) ~ (1 - q) e ^ D ^((x c - x)(l - q)^) 




(12) 




n 




(13) 
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Proposition 2. Assume that, for a class of self- avoiding polygons without punctures, the 
area moment coefficients p>m have the asymptotic form, for k G N 0; 

pj? ~ A k x- m m? k ~ x (m -> oo), (14) 

/or numbers > 0, x c > and exponents jk = {k — 9)/(j), where 6 and <fi are real constants 
and < 4> < 1 . Then, the area moment coefficient Pm of the polygon class with r > 1 
minimal punctures is asymptotically given by, for k G No, 

Pm (nfc) ~ A^x^m^- 1 (m -> oo), (15) 
waere = Ak+ r xl r /r\ and 7^ = 7A: +r .. 

Proof. We will derive upper and lower bounds on p m ' , which will be shown to coincide 
asymptotically. Let us call two polygons interacting if their boundary curves have non- 
empty intersection. An upper bound is obtained by allowing for interaction between all 
constituents of a punctured polygon. Let a polygon of half-perimeter m — 2r and area 
n + r be given. The number of ways of placing r squares inside & is clearly less than 
(n + r) r /r\. We thus have 

pT' k) < Pnt k) ■■= ^ J2 n ^ n + r ) r Pm-2r,n+r = ^ ^ fa ~ r) fc n> m _ 2r>n . 
n>l n>r+l 

By Bernoulli's inequality, we get for p~m the bound 

i_ J2 {n k+r -km^-^Pm^KpfrV <^ ^Pm^n. 
n>r+l n>r+l 

For every polygon of perimeter 2s and area t we have t > s — 1. Thus, for m sufficiently 
large, we can replace the lower bound of summation r + 1 by zero. In particular, the latter 
relation is for m > 3r + 2 equivalent to 

1 f (fc+r) _ i._._(A+r-l)\ ^ ~(r,fc) / 1 Jfc+r) 



^m-2r Pm-2r J — Pm — 



m— 2r ' 



The assumption Eq. (|T4"|) on the asymptotic behaviour of then implies that 

rl 

We derive a lower bound by subtracting from the upper bound an upper bound on the num- 
ber of square-square and square-boundary interactions. Clearly, square-square interactions 
are only present for r > 1. For a given polygon the number of square-square interac- 
tions of r squares is smaller than the number of interactions between two squares, where 
the remaining r — 2 squares may occur at arbitrary positions within the polygon. There 
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are five possible configurations for an interaction between two squares, yielding the upper 
bound 5(n + r)(n + r) r ~ 2 . Thus, the contribution to from square-square interactions 
is bounded from above by 

^n fc 5(n + r)(n + r) r - 2 p m _2r,n+r = 5(r- l)!p&" 1,fc) , 

n>l 

which is asymptotically negligible compared to p"m k \ Similarly, the number of configura- 
tions arising from square-boundary interactions is bounded from above by Y^j=i ^ ( m ~ 
2rY(n + r) r_J . This bound is obtained by estimating the number of configurations of j 
squares at the boundary by 4 J (m — 2r)- J , the factor 4 arising from edge and vertex inter- 
actions, the factor (n + r) r ~i accounting for arbitrary positions of the remaining (r — j) 
squares. We thus get an upper bound 

r r 

4 j (m - 2r) J (r - ])\p~t~ j ' k) ~ Yl ^^- 2 \m - 2r)^ +r - J ' ) 

3=1 3=1 

~ Ax 2 ;- 2 mp^ +r - l) (m->oo). 

By assumption, the latter bound is asymptotically negligible compared to p£i' k \ Thus, the 
lower bound is asymptotically equal to the upper bound, which yields the assertion of the 
proposition. □ 

Remarks, i) Proposition [2] expresses the asymptotic behaviour of the area moment coef- 
ficients of minimally punctured polygons in terms of those of polygons without punctures. 
The assumption Eq. ( 1T41) on the growth of the area moment coefficients of the model with- 
out punctures is satisfied for the usual polygon models [35]. The asymptotic behaviour of 
some models satisfying = 1, to which Proposition [2] does not apply, has been studied in 



ii) As discussed in the previous subsection, the amplitudes A k are related to the amplitudes 
fk of Eq. (jSJ) by Eq. if 7 fc G ( — 1, oo) \ {0}. For staircase polygons, where 9 = 1/3 and 
<p = 2/3, we have explicit expressions for the amplitudes A^. More generally, it has been 
shown [T3], (33J EH that, for classes of polygon models whose generating function satisfies 
a g-functional equation with a square root as the dominant singularity of their perimeter 
generating function, we have fk = Ckfi fl~ k , where the numbers c& are, for k > 1, given by 



k 

1 , 

Ik-lCk 



1 k 

1 + 7^ c fe -jQ = 0, c = 1. (16) 

1=0 

The critical point x c as well as /o and fx are model dependent constants. For staircase 
polygons we have x c = 1/4, f = — 1 and f\ = —1/64. 

Hi) Rooted self-avoiding polygons are conjectured to also have the exponents 9=1/3 and 
= 2/3. In this case the asymptotic form Eq. (1T4"1) and the form of the amplitudes Ak, 
given in Eqs. (|3j) and (fT6|) . has been tested for k < 10 and shown to hold for to a high degree 
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of numerical accuracy [TJ]. Here x c = 0.14368062927(2) is the radius of convergence of the 
(rooted) SAP half-perimeter generating function, f = —0.929607(1) and fi = —x c /(8tt) 
are the rooted SAP critical amplitudes as in Eq. ([6]). We conjecture that the asymptotic 
form (THl) holds for rooted SAPs for all values of k. Accepting this conjecture to be true, 
Proposition [2] gives the asymptotic behaviour for rooted self-avoiding polygons with r 
minimal punctures. By definition, unrooted SAPs have the same amplitudes A^. 
iv) The crude combinatorial estimates of interactions in the proof of Proposition [2] cannot 
be used to obtain corrections to the asymptotic behaviour. See also the discussion in the 
conclusion. 

4.2 Polygons with r punctures of bounded size 

The arguments in the above proof can be applied to obtain results for polygon models 
with a finite number of punctures of bounded size. The following theorem generalises 
Proposition [2] and serves as preparation for the next section, where the case of a finite 
number of punctures of arbitrary size is discussed. For a class of punctured self-avoiding 
polygons, consider their area moment coefficients 

(r,k,s) ._ ST^ k (r,s) 
fm • / , " fm^n ' 

n 

where pmfi denotes the number of polygons in the class of half-perimeter m and area n 
with r punctures, r e N , obeying the condition that the sum of the half-perimeter values 
of the puncturing polygons equals s. For simplicity of notation, we write p m ^ n := 

Pm := Pm k,0) and p m := pm . 

Theorem 1. Assume that, for a class of self-avoiding polygons without punctures, the area 
moment coefficients p£} have the asymptotic form, for k € No, 

pt ] ~ A k x~ m m? k - X (m -> oo), 

for numbers A^ > 0, x c > and 7^ = (k — 9)/<p, where 9 and are real constants and 
< (p < 1. Denote its half-perimeter generating function by V(x) = (^2 m>0 x m Pm)- Fix 

r > 1 und s6N such that [x s ](V(x)) r 7^ 0. Then, the area moment coefficient pm k ^ of 
the polygon class with r > 1 punctures whose half-perimeter sum equals is asymptotically 
given by, for k EN , 

Pm M ~ A^x^m^- 1 (m -> 00), (17) 
where 7^ = ^k+r and A k r ^ = ^ JL x s c [x s ](V(x)) r . 

Remarks, i) With s — 2, Theorem [T] reduces to Proposition [2J By summation, we also 
obtain the asymptotic behaviour for models with r punctures of total half-perimeter less 
or equal to s. Note that we have the formal identity 

00 

s=0 
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The above expressions are convergent for |x| < x c . If 6 > 0, the sum is also convergent in 
the limit x /* x c . 

ii) The remarks following the proof of Proposition [2] also apply to Theorem [H 

Proof. This proof is a direct extension of the proof of Proposition [2] to the case of a finite 
number of punctures of bounded size. We consider a model of punctured polygons where, 
for fixed s, the r punctures of half-perimeter Sj and area U satisfy s\ + . . . + s r — s. We give 
an asymptotic estimate for pm k ' s ^ ■ Let a polygon 8? of half-perimeter m — \s\ and of area 
n+ \ t\, where \s\ = s\ + . . . + s r and |t| = t\ + . . . + t r , be given. To obtain an upper bound 
for pm k ' s \ ignore all interactions between components of a punctured polygon. Recall that 
two polygons interact if their boundary curves have non-empty intersection. The number 
of ways of placing r punctures inside & is clearly smaller than 

(n+\t\) r /r\. 

This bound is obtained by considering the number of ways of placing the lower left corner 
of each puncture on each square plaquette inside the polygon. Note that, unlike in the 
proof of Proposition [2j this bound also counts configurations where punctures protrude 
from the boundary of . We will compensate for these over-counted configurations when 
deriving a lower bound for pm*' 8 '. We have 

<P^ k ' S) := ^E| s |= s E^E„>l™ fc (™+ \t\YPm-\ S \,n+\t\U r i=lPsi,ti 

= H E| S |= S E tj E„>[i|+l(l - W^Pm^slnU^lP^ta (18) 

where the first sum is over the variables S\, . . . , s r subject to the restriction \s\ = s, and 
the second sum is over all values of the variables t%, . . . , t r . Note that, for m fixed, all sums 
are finite. Invoking Bernoulli's inequality, we obtain the bound 

^EE En>| t|+a K +r - m PmHs| , n uup^ < p± As) 

' \s\=s U 

— H E[ s |=s Et 4 En>|t|+1 71 + Pm-\s\,nYli=lPsi,ti- 

Consider first the asymptotic behaviour of the expression 

r 

\s\=s U n>|t[+l i=l 

If m > |s| 2 + |s| + 2, then the lower bound of summation on the index n may be replaced 
by zero. This follows from the estimate £j < sf, being valid for every self-avoiding polygon 
of half-perimeter and area t,. Thus |t| < |s| 2 , and we argue that n > m — \s\ — 1 > 
|s| 2 + l>|t| + l. We thus get for m sufficiently large 

\s\=s i=l \\s\=s i=l 



[m 



oo 
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where the sum in brackets is finite. We now analyse the second term in the estimate derived 
from the Bernoulli inequality. To this end, define 

r 

\a\=s ti n >\t\+l i=l 

Using the estimate |t| < \s\ 2 , we get 

bm, s < (epSh.t i) iim ~ (e^ s| ) (™ - ~)- 

Now set 6 mjS := 6 m>s / ' {x~ m vrC^ The above estimate yields linim^oo 6 mjS = 0, since 

< < 1. 

We now derive a lower bound for by subtracting from p^'^ an upper bound on 

the contributions arising from puncture-puncture interactions and from puncture-boundary 
interactions. We will show that the lower bound coincides asymptotically with the upper 
bound, which then implies the assertion of the theorem 



r 



\\s\=s i=l 



For any polygon the number of puncture-puncture interactions between r > 1 punctures 
is smaller than the number of puncture-puncture interactions of two punctures with the 
remaining r — 2 punctures occuring at arbitrary positions in the polygon. We thus get the 
upper bound 

(ti + 4si)t 2 (n + \t\)(n + \t\) r - 2 < 6tMn + \t\y~\ 

where we used ti > s± — 1. The factor (ti + 4si)t 2 bounds the number of configurations 
of two interacting punctures, and the factor (n + |t|) r-2 arises from allowing arbitrary 
positions of the remaining r — 2 punctures. Define 

r I r 

c m ,s ■= Yl J2J2tit 2 n k (n+ \t\) r - 1 p m -\ aln+ \ t \\{ps i ,t i < s 4 ^pt-U^XiPsi 

\a\=s U n i=l \\s\=s i=l 

(r) 

where we used ti < \t\ < \s\ 2 for the last inequality. Setting c mjS := c myS /(x~ m m' y k 

we infer that linv^^ c mjS = 0. We have shown that for fixed s the puncture-puncture 

interactions are asymptotically irrelevant. 

We finally estimate the puncture-boundary interactions. This is done similarly to the 
above treatment of puncture-puncture interactions. The number of puncture-boundary 
interactions is bounded from above by 



(m - \s\) j si • . . . • Sj (n + \t\ 

3=1 



,r-3 
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where j punctures interact with the boundary, each contributing a factor 4(m— |s|)sj, and 
r —j punctures have arbitrary positions, each contributing a factor (n+ \t\). Note that the 
over-counted configurations in pyn' k ' s \ which protrude from the boundary, are compensated 
for by the above estimate. Define 

d m ,s ■= E| s |= s E tl EJ m - \s\) j n k (n+ \t\y- j Sl Sjp m -\ B \ in+ \n UH =1 p 9i ,u 

Defining d m:S := d miS /(x~ m m 7 fc we infer that linin^oo d m , s = 0. We have shown 
that for fixed s the puncture-boundary interactions are asymptotically irrelevant. This 
completes the proof. □ 



4.3 Polygons with r punctures of arbitrary size 

For a class of punctured self-avoiding polygons, consider for k e No their area moment 
coefficients 

(r,k) . = V- fc (r) 
fm ■ / j 11 ym,w 

n 

where pm,n '■= E^Lo^™'™ < 00 denotes the number of polygons in the class of half-perimeter 
m and area n with r punctures of arbitrary size, r e No- For simplicity of notation, we 
write p m ,n — Pm,n and = pm' k K In the sequel, we will use the area moment generating 
functions Vk{x) = ^pm'i™ of the model without punctures. 

Theorem 2. Assume that, for a class of self-avoiding polygons without punctures, the area 
moment coefficients p$ have the asymptotic form, for k e N 0; 

p (k) ~ A fc a; c - m m^- 1 (m -> oo) 

for numbers Ak > 0, x c > and '-fk = {k — 0)/<f>, where < < 1. Le£ Vk{x) = J^P^x" 1 
denote the k th area moment generating function. 

Then, the area moment coefficient pm^ of the polygon class with r > 1 punctures of 
arbitrary size is, for k G N , bounded from above by 



P 



(r>Jfc) ^g^K^g 

r! 



For /imte critical perimeter generating functions, characterised by 6 > 0, p m ' is asymp- 
totically given by, for G No, 



,(r,fc) 



[x m }V k+ r(x)(V (x)Y 



fm 



A^x- m m^- 1 (m -> oo), (19) 



where the amplitudes A^ are given by 



A (r) = A k+r (Vo(x c ))^ (2Q) 
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where Vq{x c ) := \im x y Xc Vo(x) < oo is the critical amplitude of the half-perimeter generat- 
ing function. 

Remarks, i) The asymptotic form Eq. ( |T9j) is formally obtained from Theorem [1] in the 
limit of infinite puncture size, see Remark i) after Theorem [U This observation is also 
the main ingredient of the following proof, by noting that the upper bound has the same 
asymptotic behaviour. 

ii) For staircase polygons, where 9 = 1/3 and = 2/3, the assumptions of Theorem [2] are 
satisfied. For self-avoiding polygons, we have the numerically very well established values 
9 = 1 and <fi = 2/3, which we believe to describe the asymptotic behaviour of SAPs. For 
models satisfying 9 < 0, the upper bound generally does not coincide asymptotically with 

(r k) 

Pm ■ An example of failure is rectangles with a single puncture. 

Proof. We obtain as in the proof of Theorem [1] an upper bound pm ^ for the area moment 
coefficients pm k \ It is given by 

s=0 \ s \=s i=l 

Assume in the following that 9 > 0. The asymptotic behaviour of the rhs of ( |T9l) follows by 
r-fold application of Lemma 1, which is given in the appendix. Note that, for M arbitrary, 
we have by definition 

M 

(r,k) > n (r,k,s) 

rm — / j "m ' 
s=0 

where pm' k ' s ^ is the number of r-punctured polygons, whose punctures have total perimeter 
equal to s. Theorem [T] implies that the above sum is, for M sufficiently large, asymptotically 
in m, arbitrarily close to the upper bound p& k \ See also the remark following Theorem [TJ 
This yields the statement of the theorem. □ 

4.4 Limit distribution of area and scaling function conjectures 

We first discuss the implications of the previous results on the asymptotic area law of 
polygon models with punctures. By an application of Proposition [TJ Theorem [1] and 
Theorem [2] immediately yield the following result: 

Theorem 3. Assume that, for a class of self-avoiding polygons without punctures, the area 

(k) 

moment coefficients p m have the asymptotic form, for k G N 0; 

p^ ~ A k x- m m^- x (m -> oo) 

for numbers A\~ > 0, x c > and 7/- = (k — 9)/cf), where < <fi < 1. Assume further that 
the numbers A k satisfy the Carleman condition 

fc>0 
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Denote the half-perimeter generating function of the model by V{x) = (Y2 m >o x7n Pm) ■ 

i) Consider for r > 1 the corresponding model with r punctures of bounded size, whose 
half-perimeter sum equals s G N, such that [x s ](V(x)) r 7^ 0. Denote the random 
variables of area in the uniform fixed perimeter ensemble by Xm ■ Then, we have 
convergence in distribution, 



[m — » 00 



m l/<t> 

for a uniquely defined random variable X^'^ with moments 

E[(X^) k ] = ^— , 

where the numbers A^' s ^ are those of TheoremUl We also have moment convergence. 

ii) Let Xm denote the random variable of area in the fixed perimeter ensemble for the 
model with r > 1 punctures of unbounded size. If 9 > 0, then we have convergence 
in distribution, 

% — * (r) (m - 00) 

771 

for a uniquely defined random variable X^ with moments 

A (r) 

E[(X^) fc ] = 

where the numbers A£' are those of Theorem® We also have moment convergence. 

Hi) If 9 > 0, the random variables X^ and X^'^ are related by 

x s c [x s ](V(x)) r X^ = (P(x c )) r X^ s \ 

where V(x) is the half-perimeter generating function of the polygon model without 
punctures, and where V(x c ) = lim x y Xc V(x) < 00. 

□ 

Remarks, i) For a given polygon model satisfying the assumptions of Theorem [3], the 
area moments satisfy asymptotically 

E[( *f )fc] ~ D^m^ (m-00), 

for positive numbers For classes of polygon models whose generating function satisfies 
a g-functional equation with a square root as the dominant singularity of their perimeter 



20 



Amplitude 


r = 




r = 1 




r = 2 




D 2 /D\ 

D 3/ D t 


0.530518x10" 


u 


0.530143x10" 


u 


0.529356x10" 


-u 


r\ 1 noA/i /i . , 1 a- 

0.198944x10 


-0 


0.198369x10 


-0 


0.197361x10 


-0 


D 4 /Df 


0.592379x10" 


-1 


0.588127x10" 


-1 


0.581533x10" 


-1 


D 5 /D* 


0.149079x10" 


-1 


0.146994x10" 


-1 


0.144042x10- 


-1 




0.329453x10" 


-2 


0.321705x10" 


-2 


0.311511x10- 


-2 


D 7 /Dl 


0.655743x10" 


-3 


0.632288x10" 


-3 


0.603260x10- 


-3 


Ds/Df 


0.119654x10" 


-3 


0.113600x10" 


-3 


0.106501x10- 


-3 


Dg/Df 


0.202754x10" 


-4 


0.189015x10" 


-4 


0.173673x10" 


-4 


D w /Dl° 


0.322150x10- 


-5 


0.294132x10- 


-5 


0.264251x10" 


-5 



Table 1: Universal amplitude ratios for staircase polygons with r punctures. 



generating function, the amplitude ratios -DjfV D[ n are universal, i.e., independent of 

the constants fo, f\ and x c , which characterise the underlying model [T3l 1331 134"] . This 
follows from Eqs. ([5]) and ( TT61) by a straightforward calculation. The numbers are listed in 
Table [T]for small values of r. Note that the same numbers appear for punctures of bounded 
size. 

ii) For the above class of models, explicit expressions for the asymptotic behaviour of 
their moment generating functions and their probability distributions can be derived from 
the area amplitude series via inverse Laplace transform techniques. Since the resulting 
expressions are quite cumbersome, we do not give them here. For r = the corresponding 
limit distribution of area is the Airy distribution [28l [291 EH] . The extension to r > 1 is 
straightforward. As mentioned above, for r = the amplitude ratios are found to coincide 
with those of (rooted) self-avoiding polygons to a high degree of numerical accuracy [14"] . 
If the conjecture holds true that they agree exactly, then the above expressions for limit 
distributions also appear for rooted self-avoiding polygons, for all values of r. See Section 
6 for a detailed numerical analysis. 

We now discuss the relations between the area amplitude series F(z) of the polygon 
model without punctures and F^ r '(z) of the polygon model with r punctures. Since all of 
our models have an entire moment generating function, the Carleman condition is satisfied, 
and Theorem [1] and Theorem [2] yield, by a straightforward calculation, the following result. 

Theorem 4. Assume that, for a class of self-avoiding polygons, the polygon model without 
punctures has an area amplitude series, given by 

k>0 

where 7^ = (k — 6)/<p G ( — l,oo) \ {0} and < <f> < 1, and where the numbers fk 7^ 
are related to the amplitudes in Proposition U\ via Eq. (TJJ). Denote the half-perimeter 
generating function of the model by V{x) = J2 m>0 x m p m . 



,(r) 



1 k 
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i) Assume that r > 1 and s G N are given such that [x s ](V(x)) r 7^ 0. Then, the 
corresponding model of punctured polygons with r punctures of bounded size s has an 
area amplitude series, given by 

k>0 

where 7^ = (k — 8 r )/<j) r . We have 6 r = 9 — r, r = cf), and 

F^\z) = t^ <[x s mx)r J2(k)J k z-^. (21) 

k>r 

ii) If 6 > 0, the corresponding model of punctured polygons with r > 1 punctures of 
arbitrary size has an area amplitude series, given by 

k>r 

where V(x c ) = \im x y Xc V(x) < 00. 

□ 

Remarks, i) Eq. (|2T|) allows one to derive explicit expressions for the area amplitude 
series in terms of F(z). For models where 9 = 1/3 and = 2/3 such as staircase polygons, 
the area amplitude series F(z) satisfies the Riccati equation 

F{zf - AhF\z) - f*z = 0. (22) 

This can be used to show that F^^^z) (and also F^(z)) is of the form 

r+l 

F^){ z ) = Y,PkMF{z)\ 

k=0 

where Pk,r(z) are polynomials in z of degree not exceeding |~3r/2~|, and p r+ i )r (z) is not 
identically vanishing. Simple expressions for the polynomials Pk,r{z) are not apparent. 
We note, however, that such expressions appear as correction-to-scaling functions of the 
underlying polygon models without punctures [13] . 

ii) The model of rooted self-avoiding polygons has been found numerically to have the 
same type of area amplitude series as staircase polygons (with different constants fo and 
/1). Similar considerations apply to the model of unrooted self- avoiding polygons starting 
from Eq. ©. 

We finally discuss the scaling function conjectures implied by the results of the previous 
subsections. For staircase polygons, the area amplitude function satisfies the differential 
equation (|22|) . This differential equation has a unique solution T{z) analytic for ^t(z) > 0, 
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having F(z) as an asymptotic expansion at infinity, T{z) ~ F{z) as z — > oo. The function 
.F(z) is explicitly given by 



and this function coincides with the scaling function of staircase polygons Eq. (jSJ). 

In analogy to the above observation, we conjecture that the area amplitude series for 
punctured staircase polygons determine their scaling functions. Likewise, we conjecture 
that the area amplitude series for punctured rooted self-avoiding polygons determine their 
scaling functions. 

Conjecture 1. Let r > 1 and s > 2 be given. For staircase polygons and rooted self- 
avoiding polygons, the area amplitude series F ( - r ' s \z) and F <yT \z) of Theorem^ uniquely 
define functions J-^ r,s '{z) and J-v\z) analytic for ^t(z) > Zq and non-analytic at z = Zq, for 
some negative real number zq < 0. We conjecture that the functions J-^ r,s '(z) : (zq, oo) — > R 
and T^ r '{z) : (z , oo) — > R are scaling functions as in Definition® 



Remark. The above conjecture has the following implications. 

i) Staircase polygons with a single minimal puncture specialise to Eq. (TT2|) . 

ii) Up to constant factors, the scaling form of the model with r punctures is obtained as 
the r th derivative w.r.t. q of the scaling form of the model without punctures, as can be 
proved by induction. As derivatives can be interpreted combinatorially as marking, this 
reflects the fact underlying the proofs in this section that punctures may be regarded as 
being asymptotically independent, and boundary effects do not play a role asymptotically. 
Hi) Ignoring questions of analyticity, a (formal) calculation yields that the functions 
(and JF*^' 5 )) lead, for both staircase polygons and (unrooted) self- avoiding polygons, to the 
same critical exponents in the branched polymer phase as those conjectured previously 
[21| [16] . These are obtained from the singular behaviour of about the singularity of 
smallest modulus on the negative axis, i.e., at the first zero of the Airy-function on the 
negative axis, see [131 Sec 1]. The fact that V^(x,q) is obtained from V(x,q) by r-fold 
differentiation yields the result. 

5 Computer enumerations 

Here we briefly outline which algorithms were used to derive the series expansions for the 
area moments of punctured polygons. In most cases (SAPs and punctured staircase poly- 
gons) the algorithms are simple generalisations of previous algorithms already described in 
detail in other papers, referenced below. In these cases we give brief details of the length 
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Figure 3: Illustration of the transfer matrix boundary line and local updating rules. 



of the series and the amount of CPU time used. Only in the case of staircase polygons 
with minimal punctures did we write a new specific algorithm which we shall describe in 
some detail. 

The series for punctured self-avoiding polygons were calculated using a simple general- 
isation of the parallel version of the algorithm we used previously to enumerate ordinary 
SAPs [38]. In each case (SAPs with one or two minimal punctures and SAPs with one or 
two arbitrary punctures) we calculated the area moments up to k = 10 for SAPs to total 
perimeter 100. Since the smallest such SAPs have perimeter 16 and 24 this results in series 
with 42 and 38 non-zero terms, respectively. The total CPU time required was about 5000 
hours for each of the once punctured SAP problems and up to 3000 hours for the twice 
punctured problems. The bulk of these calculations were performed on the old facility of 
the Australian Partnership for Advanced Computing (APAC), which was a Compaq Server 
Cluster with ES45 nodes with 1GHz Alpha chips (this facility has since been replaced by 
an SGI Altix cluster). 

In [22] we used a very efficient algorithm to enumerate once punctured staircase poly- 
gons. The algorithm is easily generalised in order to calculate area moments which we 
have done to perimeter 718 (k = 1), 598 (k = 2) and 506 (k = 3 to 10). It is also quite 
straight-forward to generalise the algorithm to count twice punctured staircase polygons 
and in this case we obtained the series to perimeter 502 for k — 0, perimeter 450 for k = 1 
and 2 and to perimeter 302 for k = 3 to 10. It is also easy to extract data for staircase 
polygons with punctures of fixed combined perimeter. In each case we used around 1000 
CPU hours on the APAC Altix cluster which use 1.6GHz Intel Itanium 2 chips. 

Finally we describe the algorithm used to enumerate minimally punctured staircase 
polygons. The algorithm is based on so-called transfer matrix techniques. The basic idea 
is to count the number of polygons by bisecting the lattice with a boundary line. In the left 
panel of Fig.[3]we show how such a boundary (the first medium thick line) will intersect the 
polygon in several places. The first and last occupied edges intersected by the boundary line 
are the directed walks constituting the outer staircase polygon. The other occupied edges 
(if any) belong to the minimal punctures. In a calculation to maximal half-perimeter m we 
need only consider intersections with widths up to w = m/ 2. Any intersection pattern (or 
signature) can be specified by a string of occupation variables, S = {ctq, <7i, . . . cr w }, where 
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0£ = 0, 1 or 2 if edge number i is empty, an occupied outer edge or an occupied edge part 
of a minimal puncture, respectively. We could use the same symbol for all occupied edges 
but it is convenient to explicitly distinguish between the two cases. For each signature we 
keep a generating function which keeps track of the number of configurations (to the left of 
the boundary line), that is, it counts the number of possible partially completed polygons 
with a given signature. In order to count the total number of punctured polygons we move 
the boundary line to the right column by column with each column built up one vertex 
at a time. In the left panel of Fig. [3] we have also shown a typical move of the boundary 
line, which starts in the position given by the second medium thick line and where we 
add two new edges to the lattice by moving the kink in the boundary line to the position 
given by the thin lines. As we move the boundary line to a new position we calculate the 
associated generating functions (the updating rules will be given below). Formally we can 
view this transformation between signatures as a matrix multiplication (hence our use of 
the nomenclature transfer matrix algorithm) . However, as can be readily seen, the transfer 
matrix is extremely sparse and there is no reason to list it explicitly (it is given implicitly 
by the updating rules). 

We start the calculation with the initial signature {1, 1, 0, . . . , 0}, which corresponds 
to inserting the two steps of the outer walks in the lower left corner (the count of this 
configuration is 1). As the boundary line is moved it passes over a vertex and the updating 
depends on the states of the edges below and to the left of this vertex. After the move 
we 'insert' the edges to the right and above the vertex. There are four possible local 
configurations of the 'incoming' edges as illustrated in the middle panels of Fig. [3j Both 
edges are empty, one of the edges is occupied and the other edge is empty or both edges 
are occupied. 

Both edges empty: If both incoming edges are empty then both outgoing edges can be 
empty. Else we may insert two new steps which must be part of a minimal puncture 
(the outgoing edges are in state '2'). This is only possible if the vertex is in the interior 
of the polygon (there is an edge in state '1' both below and above the vertex). 

Left edge empty, bottom edge occupied: The walk occupying the incoming edge 
must be continued along an outgoing edge. If the occupied edge is part of the 
external polygon (in state '1') there are no restrictions. If the occupied edge is part 
of a minimal puncture the walk can only be continued along the edge to the right of 
the vertex (otherwise we would not get a minimal puncture). 

Left edge occupied, bottom edge empty: This is similar to the previous case except 
that an edge in state '2' must be continued along the edge above the vertex. 

Both edges occupied: If both edges are in state '2' we close the puncture and the 
new edges are empty. If the incoming edges are in state '1' we have closed the outer 
polygon and then we add the count to the running total for the generating function. 

In the last panel of Fig. [3] we show how the updating rules given above through a sequence 
of moves of the boundary line gives rise to a minimal puncture. 
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The perimeter of the completed polygon is given by the position of the boundary line 
when the polygon is closed, e.g., if we have taken k steps in the x-direction (completed k 
columns) and moved the kink I steps in the ^/-direction then the outer half-perimeter is 
k + I and the total half-perimeter in k + I + 2r. So we need only keep track of the number 
of punctures r. This is done by associating a truncated polynomial Ps(x) = ^2p r x r with 
each signature, where the coefficient p r is the number of partially completed polygons 
with r punctures of the signature S. As a new signature S' is created from S we set 
Ps'(x) = Ps'(x) + x 5 Ps(x), where 5 = 1 if an additional puncture is inserted (as in the first 
case described above) or otherwise. The extension to the calculation of area moments is 
described in |39j. 



6 Numerical analysis 

We now turn to our numerical analysis of the series for punctured polygons. In section 16.11 
we use our series to determine numerically the exact area moment generating functions for 
minimally punctured staircase polygons with up to 5 punctures and k < 10. The resulting 
exact expressions for the leading amplitudes are in complete agreement with the formula 
derived in Proposition [2j In section 16721 we extend the study to staircase polygons with one 
puncture of fixed size and two punctures of fixed combined size. Again we find the exact 
generating functions and confirm the formula for the amplitude given in Theorem [U Next, 
in section 16.31 we analyse area moments for staircase polygons with a single puncture of 
arbitrary size. Guided by results obtained from an analysis of the conjectured exact ODE 
[22] satisfied by the perimeter generating function we carry out a careful numerical analysis 
of the area moment coefficients. This allows to obtain accurate estimates for the leading 
amplitudes and we confirm the results of Theorem [2] to at least 15 significant digits. Then, 
in section 16.41 we extend our study to staircase polygons with two staircase punctures of 
arbitrary size and we again find good agreement with the exact results. Intriguingly, we 
find in all of the above cases, that the amplitude of the leading correction term is a constant 
times the corresponding amplitude with one less puncture. Finally, in section l6"31 we present 
the results of our analysis of self-avoiding polygons with one and two punctures (minimal 
as well as arbitrary). In this case the numerical evidence is not quite as convincing, but 
we do find that the numerical estimates agree with the exact formulas to at least 3-4 
significant digits. 

6.1 Staircase polygons with minimal punctures 

In [16] it was found that the half-perimeter generating function of staircase polygons with 
a single minimal puncture is: 

mnm , N 1 - 8x + 20x 2 - 16x 3 + 2x 4 1 - 6x + 10x 2 - 4x 3 

V 1 (X) = ; r ; . 23 

1 ; 2(1 -Ax) 2 V / T^4^ 1 ; 
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This result is also derivable from Eq. (Ilip . which gives a functional equation for the area- 
perimeter generating function. It is thus plausible to expect that the generating function 
of staircase polygons with r minimal punctures is of a similar form 



where A r (x) and B r (x) are polynomials and 7 r = (3r — l)/2. We find this to be correct 
for all the cases we have enumerated that is up to r = 5: 



2A 2 {x) = x 2 - 26x 3 + 228x 4 - 906x 5 + 1709x 6 - 1378x 7 + 322x 8 , 

2B 2 (x) = -x 2 + 24x 3 - 182x 4 + 586x 5 - 815x 6 + 404x 7 - 32x 8 . 
A 3 (x) = x 2 - 22 x 3 + 197 x 4 - 924 x 5 + 2545 x 6 - 5374 x 7 + 13828 x 8 

-33634x 9 + 46027x 10 - 24746x n + 612x 12 + 256x 13 + 256 x 14 , 
B 3 (x) = -x 2 + 20 x 3 - 159 x 4 + 642 x 5 - 1509 x 6 + 3176 x 7 - 9040 x 8 
+ 19254 x 9 - 18943 x 10 + 4968 x 11 + 768 x 12 + 256 x 13 . 

2Ai(x) = 2 x 2 - 60 x 3 + 809 x A - 6564 x 5 + 36321 x 6 - 146436 x 7 
+439283 x 8 - 960070 x 9 + 1485167 x 10 - 1823356 x 11 
+2728708 x 12 - 4441406 x 13 + 4054296 x 14 - 932228 x 15 
-298318 x 16 - 143360 x 17 + 16384 x 18 - 32768 x 19 , 

2 J B 4 (x) = -2 x 2 + 56 x 3 - 701 x 4 + 5266 x 5 - 26987 x 6 + 100694 x 7 
-276415 x 8 + 537888 x 9 - 727683 x 10 + 889018 x 11 
-1536634 x 12 + 2199158 x 13 - 1289388 x 14 - 47472 x 15 
+26880 x 16 + 50176 x 17 - 6144 x 18 + 8192 x 19 . 

2A 5 (x) = 2 x 2 - 76 x 3 + 1343 x 4 - 14776 x 5 + 114384 x 6 - 666240 x 7 
+3036602 x 8 - 11071408 x 9 + 32642310 x 10 - 77911156 x 11 
+ 148630330 x 12 - 220310536 x 13 + 250700412 x 14 
-250317844 x 15 + 290657417 x 16 - 309183568 x 17 
+ 150313538 x 18 + 21743832 x 19 - 15222464 x 20 + 449152 x 21 
-3828224 x 22 - 2844672 x 23 + 974848 x 24 - 819200 x 25 , 

2B 5 (x) = -2x 2 + 72x 3 - 1203 x 4 + 12506 x 5 - 91510 x 6 + 504084 x 7 
-2171612 x 8 + 7467208 x 9 - 20683474 x 10 + 46059704 x 11 
-80841764 x 12 + 107986392 x 13 - 111525400 x 14 
+114888220 x 15 - 142562573 x 16 + 122527230 x 17 
-24478856 x 18 - 17117496 x 19 - 533632 x 20 - 2988544 x 21 
-808960 x 22 + 401408 x 23 - 819200 x 24 . 

Likewise the generating functions for the fc'th area moment, (x), is also of the 
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form (12 



v □ (r) , s = A r ,k( x ) + B rtk (x)y/1 - 4x 

We find that the degrees of the polynomials A r ^{x) and B r>k {x) are less than 5r + 2 A; for 
r < 5 and k < 10. In particular we have: 



, □(!)/ \ _ l~14a+72a: 2 ~162a 3 +145a: 4 -34a: 5 +2x 6 



Pi ( x ) - (l-4x) 5 / 2 

r 2 -82:r 3 
(l-4x) 2 



l-12a:+50x 2 -823; 3 +433; 4 -4g s ^g) 



From these solutions we then calculate the exact leading amplitudes and indeed we find 
that 

A k = A k>r (x c )/T( lr+k ) = rlx?+rr{lk+r) • (27) 

in complete agreement with Eq. (JHJ). 

We have also looked at the amplitudes = B k)r (x c )/T('y r+k — |) of the correction 
terms and find, quite remarkably, that they are given simply in terms of the leading am- 
plitudes with one less puncture: 

Bt ] = -l4~ 1] - (28) 

6.2 Staircase polygons with staircase punctures of fixed size 

Next we examine the more general case of staircase polygons with punctures of fixed size. 
In the case of one puncture we thus look at staircase polygons with a staircase hole of 
half-perimeter s while in the case of two punctures we look at staircase polygons with two 
staircase holes whose half-perimeters sum to s. As in the previous section we expect the 
generating functions to be of the form 

v (r,s) ( s = A,s,k( X ) + B ri8jk (x)y/1 - Ax 

k {X) (l-4x)TH-* ' 1 ' 

We find this to be true with the degree of the polynomials less than 5r + 2(k + s). 

For once punctured polygons we calculated the generating functions for s < 25 and 
k < 10. In Theorem CD we proved that the leading amplitude A k = Ai^ki^c) / r (jk+i) = 
A k+ \x s c p s , where p s is the number of staircase polygons of half-perimeter s. This formula is 
naturally confirmed by our numerical results. Of more interest is the sub-leading amplitude 
B^' s \ We find that the result for minimally punctured polygons generalises to this case 
and -E> j [ 1 ' s ' ) = — bi jS A k . We also find that the integer sequence d s = 8 s_1 6i jS = 1, 5, 29, 182, . . . 
is given by the recurrence: 

8s 2 d s + (s + 3) (7a + 10)4+i - (s + 3)(s + 2)4+2 = 0,d 1 = 0, d 2 = 1. 
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We note that d s grows like 8 s so that &i )S grows no faster than a polynomial in s. 

For twice punctured polygons we calculated the generating functions for s < 10 and 
k < 10. As a consequence of theorem [1] the leading amplitude is given by A^'^ = 

A 2 ,s,k{ x c)/^{lk+2) = A k+2 x s c YX=2Ps-tPt- And again we find that B ( k 2 ' s) = -b 2 , s A k+1 , 
though in this case we haven't found a recurrence for the integer sequence 2 x 8 s ~ 2 b 2s = 
1,9,69,510,.... 



6.3 Staircase polygons with a single staircase puncture 

We now turn to the analysis of the area moments of staircase polygons with a single stair- 
case puncture of arbitrary size (1-punctured staircase polygons for short). In a recent paper 
[22] we reported on work which led to an exact Fuchsian linear differential equation of order 
8 apparently satisfied by the half-perimeter generating function, V^\x) = J2 m >oP 



X'' 



for 1-punctured staircase polygons (that is V^(x) is one of the solutions of the ODE, 
expanded around the origin). Our analysis of the ODE showed that the dominant singular 
behaviour is 

V®{x) ~ {1 MX l x) + B{x) + °^f - AX) + D(x)(l + 4*) 13 / 2 . (30) 

The functions A(x)-D(x) are regular in the disc \x\ < 1/4. In addition there were a pair 
of singularities on the imaginary axis at x = ±i/2, and at the roots of 1 + x + 7x 2 . Note 
that the absolute value of these singularities exceeds 1/4 and so their contributions to the 
asymptotic growth of the series coefficients are exponentially suppressed. 

We expect that the area moment generating functions, Vj^\x) } should have a similar 
critical behaviour to that of fl30|) . Indeed our analysis using differential approximants [4*0] 
revealed that at x = x c = 1/4 there is a triple root with exponents —jk+i and —7^+1 + 1/2 
(twice) which is indicative that the behaviour is 

w A{x) + \B(x) + C(x) log(l - 4x)]VT=te 

>k W~ (1 _ 4x)3fc /2 + l • W 

However, the behaviour at the singularity x = X- = —1/4 is a little more complicated. 
For the first area moment we find that there is a double root with exponents 5 and 13/2, 
while for the second area moment we find a triple root at X- with exponents 7/2, 5 and 
13/2. For higher moments the behaviour is consistent with a triple root with exponents 
(13 — 3fc)/2, (10— 3k)/ 2 and (7—3k)/2. That is, the value of the leading exponent decreases 
by 3/c/2 and there is in addition a non-analytic correction to scaling with exponent 3/2. 
The upshot of this analysis is that the asymptotic behaviour of the coefficients of V^P (x) 
should be given by 



\x m 



Vi 1] {x) ~A m ($> 3fc/2 ~ J (^ + 7^ & +c i log(^)])+(-l) m ^^™ (_15+3 ^ )/2 j , (32) 
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where we have ignored the contributions from singularities in the complex plane with 
absolute values exceeding 1/4. Estimates for the amplitudes were obtained by fitting the 
coefficients pm^ = [x m }vj?\x) to the form given above using an increasing number of 
amplitudes. 'Experimentally' we find we need about the same total number of terms at x c 
and x_. So in the fits we used the terms with amplitudes a*, bi and Cj, i — 0, . . . , K and di, 
i = 0, . . . , 3K. Going only to K with the di amplitudes results in much poorer convergence 
and going beyond 3K leads to no improvement. So we use the 6K + 4 terms pm' k ^ with m 
ranging from M to M — 6K — 3 and solve the resulting system of 6K + 4 linear equations. 




1/M 1/M 

Figure 4: The relative precision of the estimates for the leading amplitude, \ao — \ / 
against 1/M for the first (top left panel), second (top right), fifth (bottom left) and tenth 
(bottom right) area moments. 

We compare the amplitude estimates to the predictions in Section 14.31 and we find that 
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the estimated leading amplitude a is given by the exact formula 



.4 



(r) 



l) k+r (k + r)\x r J k+r 
r!x? + T( 7fe+r .) ' 



(33) 



which agrees with Eq. (T2U|) since the critical half-perimeter generating, see Eq. (jTJ), for 
staircase polygons is V(x c ) — 1/4 = x c . In [22] we studied the normalised coefficients 
r m — P^m + s/^ m - Using the recurrence relations for pm (derived from the ODE) it is easy 
and fast to generate many more terms r m . We generated the first 100000 terms and saved 
them as floats with 500 digit accuracy. We found (to better than 100 digits) that the leading 
amplitude of the normalised series was So = 1024. Going back to the normalisation used 
in this paper we find ao = 1024/4 8 = 1/64 in agreement with formula ( 1331 . In figure H] we 
plot the relative precision of the estimates for the amplitude, \a — A^\/A^ against 1/M 
for some of the area moments. Recall that for k = 1 we have a series of 352 non-zero terms, 
for k = 2 we have 292 terms and for k = 3-10 we have 246 terms. In all cases (including 
for area moments not shown) the relative precision of the estimate a is better than 10~ 15 
for K = 10. For the first area moment, where we have a longer series, the precision is 
even more impressive, being better than 10~ 20 . So in all cases we can confirm the exact 
prediction ( 1331) for the amplitude to at least 15 significant digits. 
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Figure 5: The relative precision of the estimates for the amplitude, |c — Cjpl/Cjp, against 
1/M for the zeroth (leftmost panel) and tenth (rightmost panel) area moments. 



In the previous section we saw that for minimally punctured staircase polygons the 
amplitude of the correction term is just a constant times the leading amplitudes with one 
less puncture, see Eq. (1281) . It is thus natural to ask if something similar happens in the 
more general case. And indeed we find numerically that the amplitude of the dominant 
correction in fl3Tl) (the one proportional to the log-term) giving rise to the Cj-terms in ( 132]) 
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are 

Cf > = ~^A k . (34) 

In figure [5] we have plotted the relative error between the estimate of Cq and the predicted 
value Cjp for k = and 10. The estimates for the other moments are very similar with 
the accuracy of the agreement diminishing with higher moments. So in all cases fl34|) has 
been confirmed to better than 10 digit accuracy. 

For the sub-dominant correction term the amplitude B% , as approximated by the term 
bo in (132|) . is not simply related to and in fact it even changes sign as k is increased. 



6.4 Staircase polygons with two staircase punctures 

In this section we report on our analysis of the series for staircase polygons with two 
punctures of arbitrary size. Our first task is to work out the singularity structure of the 
perimeter generating function V <y2 \x) (for which we have a series with 240 non-zero terms). 
From the general considerations of Section [2] we expect a singularity at x = x c = 1/4 with 
exponent —5/2, but given the quite complicated singularity structure of V^'(x), as detailed 
in Eq. ( 1301) and below, we would expect similar complications for V^ 2 '(x). We analysed 
the series for V^ 2 \x) using differential approximants |lQj. This analysis revealed that there 
is a triple root at x = x c = 1/4 and the exponents had values —2.499(1), —2.070(5) and 
—1.78(1). So despite having a series of 240 terms it is still very difficult to pin down 
the exponents accurately. Given the values quoted above two possible scenarios present 
themselves. Either the exponents have the exact values —5/2, —2 and —2 or they have 
the exact values —5/2, —2 and —7/4. The behaviour of V^(x) would support the first of 
these scenarios and below we shall present evidence from the analysis of the asymptotic 
form of the coefficients which strongly supports this behaviour. We also find a double root 
singularity at i = i_ = - 1/4 with exponent estimates consistent with the exact values 
5 and 11/2. In addition there are several conjugate pairs of singularities in the complex 
plane. The most important of these are at x = (—1 ± iV3)/8, which has magnitude 1/4 
and thus lies equidistant from the origin to x c . So unlike the situation for once punctured 
staircase polygons we cannot ignore the complex singularities. The exponent estimate at 
this singularity is consistent with the exact value 33/2. The singularities at x = (— 1 ± 
«V / 3)/8 are the roots of the polynomial 1 + Ax + lQx 2 which would indicate that the 
generating function contains a term ~ E(x)(l + Ax + 16x 2 ) 33//2 . Finally, we find some 
singularities with magnitude greater than 1/4. There are singularities at x — (1 ± zy / 3)/6 
(which has magnitude 1/3) with an exponent 33/2 (we note that these are the roots of 
(1 — 3x + 9x 2 )) and at x — ±i/2 (magnitude 1/2) with an exponent consistent with the 
value 5. We also find weak evidence that V^(x), just as V^'(x), has a pair of singularities 
at the roots of 1 + x + 7x 2 . 

As noted above we need to include the contribution from a conjugate pair of complex 
singularities to the asymptotic form of the coefficients. In general this is not as straight- 
forward a task as dealing with singularities on the real axis. In j4T] we examined the 
generating function of self-avoiding walks on the honeycomb lattice which has a pair of 
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singularities on the imaginary axis at x — ±i/r, arising from a term of the form H{x){l + 

t 2 x 2 )~ v . This typically produces coefficients which change sign according to a + H 

pattern. This can be accommodated by including terms of the form 

~ r™™"- 1 ^(-l) L(m+J ' )/2j V™' 

j>0 

in fitting to the asymptotic form of the coefficients. 

Note that the analysis in [41] clearly demonstrated that, as is done above, one has 
to shift the sign-pattern by j when terms proportional to \jm? are considered. Terms 
arising from other complex conjugate pairs of singularities can give rise to much more 
complicated sign-patterns. In order to handle such cases we simply form the Taylor ex- 
pansion of the simplest possible term arising from the singularity and take the sign of the 
appropriate coefficient. Specifically in order to include terms proportional to l/m- 7 , when 
looking at the coefficients [x m ]V^(x), we take the sign to be the sign of the coefficient 
of x m+ ^ in the Taylor expansion of the function (1 + Ax + 16x 2 ) 33 / 2 . We use the notation 
Sign ([x m4t? ](l+4x + 16a; 2 ) 33 / 2 ) for this operation. 

The singularity structure revealed above leads us to fit the coefficients to a form, which 
is appropriate if the first scenario (exponents at x = x c = 1/4 are —5/2, —2 and —2) is 
correct 

[x m }V^ 2 \x) = 4 m (£f =0 ™ 3/2 - J [a 3 + ^ [b 3 +c 3 log(m)]] + (-l) m E?fo d i™" 6 " J ' /2 

+ Eifo e j Si g n (k"^](l+4x+16x 2 ) 33 / 2 ) m" 35 / 2+ ^ , (35) 

where we have ignored the singularities with magnitude exceeding 1/4. We also examine 
the alternative form appropriate if the second scenario (exponents at x = x c = 1/4 are 
—5/2, —2 and —7/4) is correct 



[x m }V^\x)= 4™(£f =0 m 3/2 - J + i + +(-l) m E-fo^ 

+ E'=o e i Si g n ([x"^'](l+4x+16x 2 ) 33 / 2 ) m" 35 / 2+ ^ . (36) 

In figure [6] we have plotted the resulting estimates for the leading amplitude ao, which 

(2") ^ 

we expect are given by the exact value A = m72y ^ , an d the amplitude cq of the sub- 
dominant terms from the two alternative asymptotic forms. First we focus on the sub- 
dominant terms. In the top right panel we show the estimates when fitting to the form (|35|) 
where we have a sub-dominant term proportional to mlog(m) while the bottom right panel 
shows the estimates obtained when fitting to the form (1361) where the sub-dominant term is 
m 3 / 4 (in both cases the dominant term is proportional to m 3 / 2 with a second sub-dominant 
term proportional to m). In the bottom right panel we note that the estimates for the 
amplitude Cq of the term m 3//4 seems to diverge. As K is increased rather than settle down 
we find that the slope of the curves of the estimates plotted vs. 1/M increases which is 
the opposite of what we would expect if we are fitting to the correct asymptotic form. We 
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Figure 6: Estimates of the leading amplitude ao and the amplitude Co of the sub-dominant 
term. The top panels show the results from fitting to the form (I35I) while the bottom 
panels are results from fitting to the form (136|) . The straight line is the exact value of this 
amplitude A y 
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take this as firm numerical evidence that (I36I) is incorrect. This could also explain why the 
corresponding estimates of a (bottom left panel) don't appear to converge to the predicted 
exact value. In contrast the estimates for Co of the term m log(m) (top right panel) from 
the form (|35l) do seem to converge and the slopes of the estimates plotted vs. 1/M decrease 
with K. In this case the estimates for a (top left panel) clearly can be extrapolated to a 
value consistent with the predicted exact value. Note further that the top left panel has 
a resolution along the ordinate which is a factor 5 higher than in lower right panel so the 
estimates when fitting to the form (135]) are much more tightly converged. The conclusion 
is that the numerical evidence clearly favours the asymptotic form ( 1351) and we believe this 
to be (if not entirely correct at least a very good approximation to) the true asymptotic 
form of the coefficients of the generating function V^ 2 '(x) for twice punctured staircase 
polygons. 
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Figure 7: Estimates of the leading amplitude A k of various area moments obtained by 
fitting to the form fl371) . The straight line is the exact value of this amplitude. 
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Now that we have settled the question of the correct singularity structure of V <y2 \ 



x 



we turn our attention to the analysis of the area moment generating function V k 



x) 



As 



for once punctured staircase polygons we find that the leading exponent at all singularities 
decreases by 3k/ 2. So in order to estimate the leading amplitude ao we fit to the form 
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(37) 



where for simplicity our notation suppresses the /c-dependence of the amplitudes. Recall 
that for k = 1 and 2 we have series with 214 terms and for k = 3-10 we have 140 terms. We 
compare the amplitude estimates to the predictions of Section HI The leading amplitude 



A k is given by the exact formula (1331) . In figure [7| we show the estimates for the 
leading amplitudes for area moments with k = 1, 2, 5, and 10. In all cases the amplitudes 
estimates appears to converge to the predicted exact value and agreement is found to at 
least 3 significant digits. 
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Figure 8: The ratio cq/A^ for the zeroth area moment (leftmost panel) and several different 
area moments (rightmost panel). 

(2) 

Finally we turn our attention to the amplitude C k of the dominant correction term. 
In the leftmost panel of figure M we have plotted the ratio between the estimated amplitude 
Co and the predicted value of A k for k = using several cut-offs. In the rightmost panel 
we show the same ratio but for several different moments using the cut-off K = 6. These 
plots are consistent with oc A k with the constant of proportionality being —0.212(2). 



6.5 Punctured self-avoiding polygons 

Before proceeding to the estimation of the amplitudes we briefly have a look at the critical 
behaviour of the area moment generating functions for punctured self-avoiding polygons. In 
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Degrees 




Once 


punctured k = 




[8,9,9,10] 


-0.002022 


- 0.023965z 


-0.002022 + 0.023965i 


0.518101 


[8,9,9,11] 


-0.003847 + 0.034768i 


-0.003847 - 0.034768z 


0.528973 


[8,9,9,12] 


-0.004553 


- 0.038104i 


-0.004553 + 0.038104i 


0.533388 


[8,9,10,10] 


-0.001830 


- 0.023143i 


-0.001830 + 0.023143i 


0.515916 


[8,9,10,11] 


-0.044131 




0.060194 


0.459558 


[8,9,11,10] 


-0.006365 


- 0.045090z 


-0.006365 + 0.045090i 


0.547479 


[8,10,9,10] 


-0.007054 


- 0.047578i 


-0.007054 + 0.047578i 


0.552445 


[8,10,9,11] 


-0.009398 


- 0.055135i 


-0.009398 + 0.055135i 


0.570859 


[8, 10, 10, 10] 


-0.012495 


- 0.063935i 


-0.012495 + 0.063935z 


0.595957 



Table 2: Biased estimates of the critical exponents of once punctured SAPs. 

[YE] we analysed the behaviour of (x) and found the critical behaviour to be consistent 
with Vi r \x) ~ A(x) + B{x)(x - x c )~^+- + C(x)(x - x c )-^+ r+1 / 2 , where ^ = 3j/2 - 3/2. 
From previous work [321 [38] we have very precise estimates for x c , which is indistinguishable 
from the positive root of the polynomial 581a; 2 + 7x — 13, that is, x c = 0.1436806292 . . .. 
Using this value for x c we form a K th -order biased differential approximant (DA) to V^\x) 
by matching the coefficients in the polynomials Qi(x) of degree iVj so that (one) of the 
formal solutions to the homogeneous differential equation 

K , 

^2(x-x e ) i Q i (x)(x—) i P{x) = 

i=0 

agrees with the first M = Ylii-^i + 1) series coefficients of T"\. (x). We are thus 'forcing' 
the DAs to have regular singular points at the origin and x c . The critical exponents 7j 
(j = 1 . . . K) are estimated from the indicial equation at x c (note that due to the biasing, 
x c is root of order K). If the true singular behaviour at x c implies a root of degree less 
than K we expect that the 'true' exponents will be quite well estimated and show little 
scatter while the 'surplus' exponents will show a lot of random scatter. In the following 
we always use K = 3 and denote the degrees of the polynomials Qi as [N 3 , N 2 ,Ni, N ]. 

First we look at the perimeter generating function for once punctured SAPs (for which 
we have series with 42 terms). In this case we have 71 = and we thus expect a logarithmic 
singularity at x c with a square- root correction term as argued in [16]. In Table [2] we list 
some exponent estimates obtained from the biased DAs. The exponent estimates are indeed 
consistent with the exact values 0,0, 1/2, which confirms the expected behaviour. 

In Table [3] we list some exponent estimates for the 2nd and 5th area moments of 
once punctured SAPs. The exponents support the expectations for the leading and sub- 
dominant exponent. The third exponent is of no significance-the nature of the differential 
approximant forces a third exponent to have some value, but its lack of convergence suggests 
it is not, in fact present. In Table H] we list exponent estimates for the 0th and 2nd area 
moments of twice punctured SAPs (for which we have series with 38 terms). Similar 
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Degrees 


Once punctured k 


= 2 


Once punctured k 


= 5 


[8,9,9, 10] 


-3.000291 


-2.493921 


-1.154513 


-7.500838 


-6.976201 


-5.107975 


[8,9,9,11] 


-3.000158 


-2.496645 


-1.258606 


-7.503030 


-6.927432 


-2.422194 


[8,9,9,12] 


-3.012422 


-2.421400 


-3.323742 


-7.500231 


-6.992184 


-5.461825 


[8, 9, 10, 10] 


-2.999849 


-2.502858 


-1.447159 


-7.517118 


-6.776728 


-8.010389 


[8,9, 10,11] 


-3.000452 


-2.491160 


-1.022360 


-7.500700 


-6.980331 


-5.251153 


[8,9,11,10] 


-3.000143 


-2.497148 


-1.307314 


-7.500366 


-6.988804 


-5.406296 


[8,10,9,10] 


-2.999481 


-2.511762 


-1.668527 


-7.505629 


-6.879658 


-14.20958 


[8,10,9,11] 


-2.999906 


-2.501910 


-1.441981 


-7.500535 


-6.984576 


-5.361394 


[8, 10, 10, 10] 


-2.999647 


-2.507601 


-1.578352 


-7.500946 


-6.974159 


-5.059683 



Table 3: Biased estimates of the critical exponents for the 2nd and 5th area moments of 
once punctured SAPs. 



Degrees 


Twice punctured k 


= 


Twice punctured k 


= 2 


[7,8,8,9] 


-1.504250 


-0.968766 


6.594627 


-4.499111 


-4.006801 


-2.544930 


[7,8,8,10] 


-1.504232 


-0.968829 


5.029521 


-4.495871 


-4.039199 


-3.083192 


[7,8,8,11] 


-1.504209 


-0.968910 


-60.90721 


-4.494258 


-4.048789 


-2.701322 


[7,8,9,9] 


-1.504210 


-0.968909 


12.03204 


-4.494568 


-4.047486 


-2.836670 


[7,8,9,10] 


-1.504197 


-0.968951 


32.28205 


-4.493872 


-4.051402 


-2.672813 


[7,8,10,9] 


-1.504205 


-0.968923 


22.27646 


-4.508207 


-3.915977 


-1.950665 


[7,9,8,9] 


-1.504178 


-0.969017 


18.18298 


-4.495340 


-4.040712 


-2.812621 


[7,9,8,10] 


-1.504177 


-0.969020 


18.81517 


-4.495136 


-4.041237 


-2.679530 


[7,9,9,9] 


-1.504176 


-0.969022 


18.58757 


-4.498396 


-4.012876 


-2.698591 



Table 4: Biased estimates of the critical exponents for the 0th and 2nd area moments of 
twice punctured SAPs. 
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comments apply to the three columns of exponent estimates as were made for the 2nd area 
moments. In all cases we get a clear confirmation of the critical behaviour observed in [16] . 

Next we turn our attention to estimates for the critical amplitudes. Proposition [2] and 
Theorem [2] tells us that the critical amplitude of the kth area moment of self-avoiding poly- 
gons with r (minimal or arbitrary) punctures is proportional to the critical amplitude of the 
(k + r)th area moment of unpunctured SAPs (the theorems also give the constants of pro- 
portionality). In order to test these predictions numerically we analyse in this section data 
for SAPs with one and two punctures. In all cases we look at the ratio r m = pm / Pm 
which should approach the relevant constant of proportionality. Given the critical be- 
haviour outlined above we expect that these amplitude ratios can be approximated quite 
well by the asymptotic form 

3=0 

So as in the analysis of punctured staircase polygons we obtain estimates for the leading 
amplitude a , by fitting to the above form truncated at some level K, and we then plot 
these estimates against 1/M. 



6.5.1 Minimal punctures 

According to Proposition [2] the amplitude of the k-th area moment of SAPs with r minimal 
punctures is 

4 r) = A k+r x 2 ;/r\. (38) 

We first analyse the area moments of SAPs with a single minimal puncture. In the left 
panel of figure [9] we have plotted estimates for amplitude ratio Aq/A\ with 7 < K < 10. 
The prediction for this ratio is x 2 c) which is plotted as a straight line. The estimates obtained 
with K = 9 and 10 are indistinguishable from the predicted value at the resolution of the 
plot. Note that the 'curvature' of the plotted values decreases as K increases. We take this 
to be a very clear indication that the ratio r m is very well approximated by the assumed 
asymptotic form. In the right panel we plot the estimates for the amplitude ratio A^' /A). + i 
for < k < 9 using the cut-off value K = 10. For small values of k these plots give firm 
numerical evidence for the correctness of Proposition [2J For higher values of k (8 and 9 
in particular) the evidence is not quite as firm though nothing in the plot would suggest 
a discrepancy with the predicted value. We again emphasise that in this case we have 
relatively short series of only 42 terms. 

Next we analyse the area moments of SAPs with two minimal punctures. The left panel 
of figure [TU] shows estimates for amplitude ratio A^ /A 2 with 7 < K < 10. The prediction 
for the ratio, is plotted as a straight line. Again we find that the estimates for K = 10 
are indistinguishable from the predicted value, though in this case the relative resolution 
is coarser than in the previous plot. In the right panel we plot estimates of the amplitude 
ratio A k /Ak+2 for < k < 8 using the cut-off value K = 10. Again we find that our 
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Figure 9: Estimates of the amplitude ratio /A k+ i for self-avoiding polygons with one 
minimal puncture. 

numerical analysis confirms the prediction to a high degree of confidence. Recall that in 
this case we have series with only 38 terms. 

6.5.2 Arbitrary punctures 

According to Theorem [2] the amplitude of the k-th area moment of SAPs with r arbitrary 
punctures is 

4 r) =A k+r (V(x c )) r /r\, (39) 

where V(x c ) is the critical amplitude of the half-perimeter generating function. 

The first step in our analysis is to obtain an accurate estimate of V(x c ). In [42J we 
obtained the rather imprecise estimate V(x c ) w 0.036 by evaluating Pade approximants to 
the generating function. Here we shall estimate V(x c ) directly from the perimeter data. 
We first tried the form 

M 

P™ x ™ ~ H*c) + ajM 1 ' 2 + a 2 /M + ■■■ 

m=0 

Using the first ten terms in this asymptotic expansion we found (with M = 55) that 
V(x c ) « 0.0362642, but a x w 4.28 x 10~ 9 and a 2 w -1.42" 7 , while a 3 ta -0.066. We are 
therefore quite confident that a\ = a 2 = 0. Upon further analysis we found convincing 
evidence that the correct asymptotic form in fact is 

M 

P™ x ™ ~ Hxc) + h/M 3 / 2 + b 2 /M 5 ' 2 + ■■■ 

m=0 
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Figure 10: Estimates of the amplitude ratio A\, /A k+2 for self- avoiding polygons with two 
minimal punctures. 



M_ K = 8 

45 0.036264215181387 

46 0.036264215181343 

47 0.036264215181305 

48 0.036264215181271 

49 0.036264215181240 

50 0.036264215181214 

51 0.036264215181190 

52 0.036264215181168 

53 0.036264215181149 

54 0.036264215181131 

55 0.036264215181116 



K = 10 
0.036264215181095 
0.036264215181088 
0.036264215181073 
0.036264215181060 
0.036264215181048 
0.036264215181038 
0.036264215181029 
0.036264215181021 
0.036264215181013 
0.036264215181006 
0.036264215181000 



K = 12 
0.036264215180475 
0.036264215181354 
0.036264215180915 
0.036264215180994 
0.036264215180970 
0.036264215181001 
0.036264215180972 
0.036264215180972 
0.036264215180969 
0.036264215180967 
0.036264215180962 



Table 5: Estimates of the critical SAP amplitude V(x c ). 
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Indeed, the observation that a\ and a 2 vanish follows from the assumed asymptotic be- 
haviour of p m . In Table [5] we have listed estimates for V(x c ) obtained using various values 
of M and the cut-off K in the asymptotic form. From this we confidently estimate that 
V{x c ) = 0.0362642151808(2). 
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Figure 11: Estimates of the amplitude ratio A k /A^+i for self-avoiding polygons with one 
puncture of arbitrary size. 

In figure [TP we have plotted estimates of the amplitude ratio /Ak+i for self-avoiding 
polygons with one puncture of arbitrary size. In the leftmost panel we look at the ratio 
Aq 1 /Ai using different cut-offs K. The rightmost panel shows the ratio Ap /A^ + i for area 
moments up to k = 9 using the cut-off K = 10. The straight line corresponds to the 
expected value V(x c ), using the estimate for this quantity obtained above. The estimates 
in the leftmost panel show some variation when plotted against 1/M, but in the limit 
M — > 00 the estimates appear to converge to the expected value (if the trend holds). As 
for minimally punctured SAPs we see an ever closer agreement as K is increased, again 
indicating that the assumed asymptotic form is reasonable. Obviously, as seen in the right- 
most panel, the estimates for high area moments are not as close to the expected value. 
However, given the trend in these estimates we are confident in stating that our numer- 
ical analysis is consistent with the results of Theorem [2J The agreement is particularly 
impressive bearing in mind that we analyse series with just 42 terms. 

Finally in figure [T2l we have plotted our estimates of the amplitude ratio A^ /A k+2 for 
self-avoiding polygons with two punctures of arbitrary size. The straight line corresponds 
to the expected value V(x c ) 2 /2. Again all estimates are consistent with the results of 
Theorem [21 though the numerical agreement is less convincing, but then again the series 
have only 38 terms. 
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Figure 12: Estimates of the amplitude ratio A k /Ak +2 for self-avoiding polygons with two 
punctures of arbitrary size. 

7 Conclusion 

We have rigorously analysed the effect of punctures on the area law of polygon models. 
In particular we obtained expressions for the leading amplitudes of the area moments for 
punctured polygons in terms of the amplitudes for unpunctured polygons (see Theorems [T] 
and [2]). For staircase polygons this led to exact formulas for the amplitudes. For self- 
avoiding polygons the formulas rely on an assumption about the asymptotic behaviour 
of SAPs without punctures. They contain constants not known exactly but estimated 
numerically to a very high degree of accuracy. Our analysis also led to conjectures about 
the scaling behaviour of these models (see Conjectured]). A proof of these conjectures is an 
open and difficult problem. Further numerical support of these conjectures might follow 
from an analysis of critical perimeter moments along the lines of |15j . 

The expressions for the amplitudes were thoroughly checked numerically. For staircase 
polygons with up to 5 minimal punctures and staircase polygons with one or two punctures 
of fixed size we used our series expansions to find the exact generating functions for area 
moments to order 10. Naturally, in all these cases the leading amplitude agrees with 
the proved formulas. Interestingly we find that the amplitude of the correction term is 
proportional to the corresponding leading amplitude with one less puncture. For staircase 
polygons with one and two punctures of arbitrary size a careful asymptotic analysis of 
the series for the area moments yielded very accurate estimates for the amplitudes, again 
confirming the exact formulas. Finally, we also analysed series for self-avoiding polygons 
with one and two punctures (minimal as well as arbitrary). In this case the numerical 
evidence is not quite as convincing, but we did find that the numerical estimates agree 
with the exact formulas to at least 3-4 significant digits. 

The numerical analysis also yielded explicit expressions for correction terms, see Eqs. ( 1281) . 
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(|34|) and subsection 16. 21 These correspond to contributions from puncture-boundary in- 
teractions and from puncture-puncture interactions. It would be interesting, but difficult 
to give a combinatorial proof of these results. The difficulty of any combinatorial proof 
becomes clearer when considering the recent closed form solution of a closely related model 
of punctured staircase polygons. In [23] one of us (IJ) considered a model of punctured 
staircase polygons in which the internal polygon is rotated by 90 degrees with respect to 
the outer polygon. The proofs in this paper never consider such restrictions on the place- 
ment of the internal polygon, that is, the internal polygon can be placed in any way one 
pleases. The results therefore carry over unaltered to the problem of rotated punctured 
staircase polygons. Interestingly, this means that the leading asymptotic form of the coeffi- 
cients is exactly the same for both models, any differences arising only from the correction 
terms. The dominant correction term to the half-perimeter generating function for punc- 
tured staircase polygons V{x) is oc log(l — 4a;)/\/i — Ax [22]. From the exact closed form 
solution to the half-perimeter generating function for rotated punctured staircase polygons 
T~ ) 'Rot{x) it follows that the first correction term is oc (1 — 4x)~ 3,/4 [13] • These differences in- 
dicate that combinatorial arguments for a proof of sub-dominant behaviour must be quite 
subtle! 

In this paper, we discussed models of punctured polygons, where punctures are of the 
same type as the outer polygon. More generally, punctures may be built from different 
polygon classes. In that situation, Theorem 1 holds, with the obvious modification, for any 
collection of polygon models as punctures. Theorem 2 holds, with the obvious modifica- 
tion, for any collection of polygon models as punctures, if the corresponding half-perimeter 
generating functions are finite at the radius of convergence x c of the half-perimeter gener- 
ating function of the outer polygon. This includes the model considered in [23] as a special 
case. We remark that the critical half-perimeter generating function of the outer polygon 
may be infinite. 

Our analytical results are based on the observation that puncture counting can be done 
using polygon area estimates, polygon boundary contributions being asymptotically negli- 
gible. In particular, effects of self-avoidance do not influence the results. This phenomenon 
is also expected to hold in higher dimension. Consider so-called polycubes, the generalisa- 
tion of polyominoes to three dimensions. Polycubes have been enumerated by volume up 
to 18 cubes, see |44j. Three-dimensional vesicles [25] are a subclass of polycubes, having 
no interior holes. Let a class of three-dimensional vesicles be given, counted by surface 
area, with a bounded number of vesicular holes. Assume that the asymptotic behaviour of 
the volume moments is known for the model without vesicular holes. If < < 1 and the 
critical surface area generating function is finite, then our method of proof can be adapted 
to describe the volume law of three-dimensional vesicles with vesicular holes. For the full 
model of closed self-avoiding orientable surfaces of genus zero on the cubic lattice, however, 
there is numerical evidence that 0=1, see e.g. the review in [3]. 

For models with minimal punctures, the number pm in Eq. (jl~5]) also counts the 
number of polygons winding around a fixed point of the dual lattice. This problem and 
its generalisation to M > 1 mutually avoiding self-avoiding polygons has been considered 
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previously by Cardy [16]. It would be interesting to consider whether this generalisation 
can also be treated using the above methods. This involves the analysis of polygon models 
satisfying 6 = 0. If 9 < 0, interaction terms are generally not asymptotically negligible, so 
that the above analysis does not yield asymptotically exact estimates. 

A major open question is the problem of self-avoiding polygons with an unbounded 
number of punctures. Here, Theorem [2] yields an upper bound. Let Qk(x) denote the 
k th area moment generating function for a model of punctured polygons with an arbi- 
trary number of punctures. Qk{x) is a formal power series, usually with zero radius of 
convergence. The (asymptotically exact) upper bound for models with r punctures yields 
an upper estimate for the number of punctured polygons with an arbitrary number of 
punctures. It is 

where -C denotes coefficientwise majorisation. Let Qk{x) denote the k th factorial area 
moment generating function for punctured polygons with an arbitrary number of punctures. 
It can be shown that an upper bound is given by 

r=0 

where Vk(x) is the k th factorial moment generating function, and V(x,q) is the perimeter 
and area generating function of the model without punctures. In particular, the half- 
perimeter generating function Qo(x) = Qo(x) is majorised by 

Qo(x) <£V(x,l+V(x,l)). 

Another problem touched upon in this article is punctured polygons in the fixed area 
ensemble. We gave a (non-rigorous) argument for values of the critical exponent in the 
branched polymer phase from the crossover behaviour of the tentative scaling function, 
thereby confirming previous results [21 J. In that phase, boundary effects are indeed crucial, 
such that our methods of deriving limit distributions cannot be applied in this situation. 
On the other hand, area laws are expected to be of Gaussian type, as is usually the case 
away from phase transition points. The same phenomenon is expected to occur in the fixed 
perimeter ensembles for q ^ 1. 

E-mail or WWW retrieval of series 

The series for the generating functions studied in this paper can be obtained via e-mail by 
sending a request to I.Jensen@ms.unimelb.edu.au or via the world wide web on the URL 
: / /www.ms.unimelb.edu.au/ ^iwan/] by following the instructions. 



45 



Acknowledgements 



The calculations presented in this paper would not have been possible without a generous 
grant of computer time on the server cluster of the Australian Partnership for Advanced 
Computing (APAC). We also used the computational resources of the Victorian Partnership 
for Advanced Computing (VPAC). CR and AJG would like to acknowledge support by the 
German Research Council (Deutsche Forschungsgemeinschaft) within the CRC701. IJ and 
AJG gratefully acknowledge financial support from the Australian Research Council. 



Appendix 

We analyse the asymptotic growth of a Cauchy product of sequences in terms of the 
asymptotic growth of its constituting sequences. The following lemma is an extension 
of [4"Tl Thm 2] to the case of generating functions with equal radii of convergence. Its 
proof relies on Lebesgue's dominated convergence theorem [48l Thm 1.34], which states 
conditions under which an exchange of limit and sum is allowed: With n £ N and k £ No 
let real numbers a n ± be given. Assume that lim^oo a n ^ =: £ M for all k and that for 
all k there is a bound |a n ,fc| < bk uniformly in n £ N. Assume that ^2 k>Q bk £ R. Then 
lim n ^oo X]fc>o a ™. fc = Sfc>o flfc e 

Lemma 1. Let two sequences (/ n )neN an d (flOneNo °f rea ^ numbers be given, with gen- 
erating functions f(x) = Yln>ofn xU an d 9( x ) = TlmXdQnX 11 ■ Assume that both generating 
functions have the same positive finite radius of convergence x c , < x c < oo. Assume that 
the sequences (f n ) and (g n ) satisfy asymptotically 

f n ~Ax- n n<-\ g^Bx^n 5 - 1 (n -> oo), 

for nonzero numbers /I / 0, B / 0, and for real constants 7, 5 satisfying 5 < and 
7 > 5 + 1. Assume that g(x c ) := \im x y Xc g(x) 7^ 0. Then, the Cauchy product of (f n ) and 
(g n ) satisfies 

n 

^2 fn-kQh = [x n ]f(x)g(x) ~ g(x c )f n (n -»• 00). 
Proof. Let f n := for n < and define for n £ N and for k £ No 

fn—k 

an ' k := Ax^m-i 9 "' 

Below, we derive a bound on |a n( fc| uniformly in n £ N and summable in k £ No- Then, 
Lebesgue's dominated convergence theorem can be applied to a n ^- Since g(x c ) 7^ 0, this 
yields the statement of the lemma. 

Note first that the assumption on the asymptotic behaviour of /„ implies the existence 
of a constant n £ N such that for all n > n Q + k we have the estimate 

k^- 1 



n—k 



Ax- n n^- 1 



<2x k c [l- 



n 
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Fix such n . We distinguish the three cases n < k, k < n < no + k, and n > n + k. If 
n < k, we clearly have a njk = 0. For n > n + k, the above estimate yields 

<2(l-£f 

We will first consider the case 7 — 1 < 0. If n > n + k, we have 1/(1 — k/n) < 1 + k/riQ. 
This implies 

k,*| <2 2 -^k l ^x k c \g k \=:b^. 
If k < n < no + k, we estimate similarly 

< |A| -1 max{l,a;^ }max m<no {|/ m |}(no + k) l ~^x k c \g k \ 

< ni^i,,-}^^} . { (2no) ?;]»i l9tli *; ° } =: »£■> ■ 

Define b k :— b k ^ + bf^ for A; e No. Then |a nj fc| < b k uniformly in n G N, and ^ fc>0 b k < 00. 
Summability follows from 

00 

^x^ 1_7 |^| < 00, 

fc=0 

since by assumption x^/c 1_7 |(7 fc | ~ \B\k 5 ~ J as k — > 00, where 7 — 5 > 1. 

If 7 — 1 > 0, uniform estimates are obtained along the same lines. In that situation, 
the factors (1 — k/n) 1 " 1 and (n + A;) 1 " 7 may be replaced by unity, resulting in simpler 
uniform bounds involving J2 k>0 x k c \g k \ < 00. □ 



\0'n,k\ 



fn—l 



Axz n nt- 1 



9k 



\@"n,k I 



J 71 — 



Ax c n n~i- 
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